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SECTION I 


INTRODUCTION 

A. BACKGROUND 

For modern military aircraft, manageable flight at high, near-stall 
angles of attack is a key effectiveness and safety factor. The safe 
angles of attack and sideslip attainable in level and maneuvering flight 
often represent the important differences in both offensive/defensive 
capability and accident histories of otherwise equivalent aircraft. In 
the past generation of fighter aircraft, surprise loss of control depar¬ 
ture generally has been of concern as a key maneuver-limiting factor. 
The more recent generation has emphslzed design for departure resistance 
but not necessarily prevention. Many operational pilots desire the 
capability to depart the aircraft on command as a last-ditch defensive 
maneuver, since departure Is generally violent with unpredictable 
results from the attacker's viewpoint. The key element then revolves 
around departure resistance (or susceptibility), warning, severity, and 
recoverability. Accordingly, the Issue of improved utility and safety 
at high angles of attack has been re-recognlzed as an area of fruitful 
research and development. 

The recent upsurge In Interest has resulted in numerous studies 
devoted to identifying the safety-related maneuver-limiting phenomena 
Inherent in: 

• Buffet 

• Pitchup or "dig-in" 

• Roll reversal 

• Wing rock 

• Nose slice 


• Rolling departure — snap roll 




Some of these phenomena are perturbations about a steady flight 
condition (e.g., buffet, wing rock), which not only limit tracking 
capability (as shown In Reference 1) but also may serve as warning of 
approaching "departure" from controlled flight (e.g., nose slice, roll¬ 
ing departure). The more Interesting and mysterious of these phenomena 
are the latter, and considerable effort has been expended in correlating 
the angles of attack at which such behavior occurs with various aero¬ 
dynamic stability and control parameters (e.g., References 2-4). Such 
open-loop "correlations," while useful, are not so satisfying or in¬ 
structive as the more positive identification of causal factors offered 
by the methods developed In Reference 5. 

The difficulty Is in identifying causal relationships in the complex 
Interactive pilot/vehicle situation during uncoordinated flight, where 
aerodynamic cross-coupling can compound vehicle dynamic characteris¬ 
tics. The previous (Reference 5) analysis and pilot simulation showed 
the static aerodynamic cross-coupling coefficients and N£ to contri¬ 
bute to closed-loop nose-slice departure susceptibility and severity in 
the A-7 aircraft. It represented an Initial successful and promising 
new attack on an old problem. But configuration differences among air¬ 
craft types, models and even loadings have been observed to produce 
gross differences in behavior through changed aerodynamics and inertial 
characteristics. Thus It is desirable that the same methods be applied 
to additional high-performance fighter aircraft having widely differing 
high AOA handling characteristics to see if further cause-effect rela¬ 
tionships can be identified. 

Accordingly, the stated or implied goals of this program are to: 

• Identify key design parameters that limit high- 
angle-of-attack maneuverability for contemporary 
high-performance attack and fighter-type aircraft. 

• Postulate fundamental aerodynamic and control 
system design methodologies that will alleviate 
the limiting conditions. 

• Formulate handling qualities requirements for high 
AOA maneuvering flight to be incorporated in MIL- 
F-8785C, the military flying qualities specifica¬ 
tion for piloted airplanes. 
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B. TECHNICAL APPROACH 


The research encompassed four major technical areas. The first was 
devoted to development and validation of aerodynamic models for two 
fighter aircraft having significantly different high AOA maneuver- 
limiting characteristics. The two aircraft selected, the F-4J and the 
F-14A, were found to represent almost opposite extremes In their range 
of departure susceptibility and severity. Attention was then turned to 
investigation of maneuver-limiting characteristics of each airframe and 
identification of causal relationships based on quaai-linear analysis 
at symmetric and asymmetric (0 i 4 0) flight conditions. The causal 
parameters were varied to alter the dynamic characteristics of each 
airframe at high AOA. Predicted characteristics were checked utilizing 
complete six-degree-of-freedom (6 DOF) models with nonlinear 
aerodynamics, and any differences between the quasl-linear frozen point 
analytic predictions and the nonlinear model results were resolved. 

The second task Involved development of methods, criteria, and an 
associated pilot rating scale, for evaluation of handling qualities when 
approaching controllability limits. Previous analysis and simulation 
(Reference 5) has demonstrated that the Cooper-Harper handling quality 
rating scale is inapplicable for departure/recovery flight situations. 
Therefore, a new pilot rating scale is required specifically for high 
angles of attack. 

The third task, pilot simulation, was performed to evaluate the 
influence of Intentionally varied maneuver-limiting characteristics 
(wing rock, nose slice, etc.) under different normal flying situations 
such as training flights and air combat tracking tasks; to identify key 


It is recognized that these linear analysis techniques are applied 
to phenomena which may be nonlinear in nature. However, the intent of 
the study is to Investigate in an understandable manner the conditions 
which can lead to or precipitate rapid changes in motion, as opposed to 
analyzing the fully developed large angle phenomena. Therefore, piece- 
wise linearity is appropriate — but care should be taken in attempting 
to extrapolate the results. 
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flying quality parameters in terms of departure resistance, warning, 
severity, and recovery for possible Inclusion in the flying qualities 
specification; and to exercise and refine as necessary the new pilot 
rating scale. 

The fourth task Involved assessment and correlation of results of 
the simulation with other applicable results and data; and formulation 
of high AOA maneuvering flight generalized design guides and flying 
qualities criteria for Incorporation into MIL-F-8785C. 





SECTION II 


AIRCRAFT HIGH AOA DTMAHIC CHARACTERISTICS 


This section presents condensed results of an extensive analysis of 
the high AOA dynamic characteristics of the two example aircraft. The 
F-4J is treated first, followed by the F-14A. The high AOA character¬ 
istics of each aircraft are described as determined in flight test and 
operational deployment. The analytical models are briefly summarized 
and example comparisons shown to demonstrate the match obtained with the 
actual aircraft. 

Results of linear, frozen-point dynamic analysis at a aeries of 
symmetric and asymmetric (3 i 0) flight conditions are summarized. This 
analysis provides insight to potential causal factors behind undesirable 
characteristics at high AOA, showing examples of advantages and limita¬ 
tions of such analysis when applied to situations in which aerodynamics 
are obviously quite nonlinear. In particular, this section identified 
lateral-directional static aerodynamic cross-coupling as a key factor in 
determining high AOA stability and controllability characteristics. 

A. F-4J AHALYSIS 

The F-4J model of the F-4 family was selected because a previously 
validated aerodynamic model (Reference 6) was in hand. However, all 
versions of the aircraft with an unslatted leading edge wing are 
reported to exhibit quite similar high AOA flying characteristics: 
increasing AOA produces buffet of increasing severity, increasingly 
adverse yaw, wing rock, mild pitch-up, stall, and finally departure. A 
representative trimmed lift curve with onset of these handling phenomena 
is shown in Figure 1 from Reference 7. Above about 11 deg AOA, aileron 
adverse yaw becomes significant and rudder is used for roll control. 

Wing rock appears primarily as a rolling-sideslipping motion. It 
la often described by pilots as a divergent dutch roll which is easily 
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• - MODUS 

Figure 1. Basic F-4 Trimmed Lift Curve (from Ref. 7) 


aggravated or caused by pilot lateral stick inputs. The angle of attack 
for onset and the severity actually depend upon several facto.. ,, includ¬ 
ing c.g. location, roll inertia or loading (clean vs. external stores), 
gear and flap settings, etc. It is generally more pronounced at light 
weight, with low roll inertia, and at forward c.g. The oscillations 
generally require 10 sec or more of sustained high AOA to develop. 
Amplitudes can be large enough (±30 deg $ and ±10 deg p) to interfere 
with tracking performance, or to be highly uncomfortable when escorting 
considerably slower aircraft. 

The mild pitch-up is alleviated by a pitch rate damper but does 
result in some lightening of stick forces. It has been described by 
some pilots as a "dig-in" during slow turns. Possibly it is more of a 
stall warning cue than a maneuver-limiting problem. 

Departure is described variously as nose slice, roll departure, or 
both. Its onset also appears to depend upon c.g. location and aircraft 
Inertia variations. Although shown in Figure 1 to occur at considerably 
higher AOA than wing rock, departure is not necessarily preceded by wing 
rock — especially if AOA is rapidly increased. 




1. Analytic Model 


The F-4J aerodynamic model (Reference 6) was derived from three 
excellent but separate data sources (References 8-10). It was assembled 
to model as closely as possible the F-4J sub- and transonlcally for a 
moving-base simulation used to train U.S. Navy pilots In air combat 
maneuvering. It therefore had to be a very good representation of the 
complete aircraft In high AOA maneuvering up through stall. The Instruc¬ 
tor pilots pronounced the simulation "credible and realistic of actual 
F-4J handling characteristics." 

The history of development of the complete aerodynamic model and the 
data package is presented in Part III, Appendix I. In brief, the aero¬ 
dynamic coefficients represent the following flight regime: 

0 < a < +110 deg 
-30 < 8 < +30 deg 

For this analysis the look-up table data were restricted to low Mach 
(<0.4) and e single altitude (h - 15,000). Weight and inertia charac¬ 
teristics representative of partially full internal fuel tanks, an empty 
600 gal centerline tank, and missile pylons on wing stations 2 and 8 
have been assumed. 

Conventional 6 DOF equations plus kinematic terms are employed with 
the moment equations In body centerline axes, the force equations In 
wind axes, and aircraft orientation angles in standard Euler axes. The 
nonlinear equations of motion, aerodynamic force and moment equations, 
auxiliary equations, etc., are detailed in Appendix I (Part III). 

A block diagram of the basic flight control and augmentation mech¬ 
anization is presented in Figure 2. Aircraft control Is exerted through 
horizontal s'cabllator ( 5 8ta b)» rudder (6 r ), aileron (6 fl ), and spoiler 
(6 g p) surfaces. The aileron and spoiler systems are interconnected 
such that lateral deflection of tha control stick produces downward 
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Figure 2* F-4j Flight Control System 





















deflection of one aileron and upward deflection of the opposite spoiler. 
Stability augmentation Is provided in all three axes, but it is common 
practice for the pilot to turn the roll damper off during high AOA 
maneuvering. The complex artificial feel system was not modeled In this 
analysis and simulation. 

2. Model Validation 

As a part of the checkout and acceptance test, the simulation based 
on the Reference 6 aerodynamic data was flown through various offensive 
and defensive air combat maneuvers, stalls, and departures by Navy in¬ 
structor pilots, who Indicated that it adequately represented the F-4J 
handling and performance. This provided the first gross validation of 
the aerodynamic model. 

Actual flight test traces of F-4J high AOA dynamic responses are not 
available to compare against those of the mathematical models. The 
available high AOA flight traces are from an F-4E stall/post-stall 
flight test (Reference 11). However, Reference 12 indicates that all 
hard-wing models of the F-4 have approximately the same stall/departure 
characteristics. The F-4E has a nose section approximately 5 ft longer 
than the F-4J; this could cause some difference in side force due to 
asymmetric vortex shedding at high AOA (Reference 13). In addition the 
flight test vehicle was equipped with a spin chute, had a reinforced aft 
fuselage structure to handle the spin chute loads, and la presumed to 
have offsetting ballast in the nose. Thus, Its pitch and yaw Inertia 
characteristics are quite different from the F-4J, as shown in Table 1. 
This difference is assumed to have minor Influence on the sequence of 
stall/departure characteristics with increasing AOA. The larger yaw 
inertia and longer nose (greater vortex shedding) of the F-4E could be 
offsetting departure susceptibility factors. But the Inertia differ¬ 
ences should produce measurable differences In dynamic response parame¬ 
ters such as dutch roll frequency, response to control Inputs, and 
angular rates associated with departure. For example, scaling parame¬ 
ters by the ratio of inertias gives ratios of tti gp and a> d of 1.11. 
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TABLE 1. COMPARISON OF F-4E AND F-4J 
WEIGHTS AND INERTIAS 


F-4E (Reference 11) F-4J 


<£. stores 

Empty tank 

Empty tank 

Wing stores 

Pylons 1, 2, 8, 9 

Pylons 2, < 

W (lb) 

40,000 

37,000 

c.g. (X c) 

28.1 

29.3 

I x (slug-ft 2 ) 

27,500 

23,850 

Iy (8lug-ft 2 ) 

157,000 

127,400 

I 2 (alug-ft 2 ) 

180,600 

146,000 

I xz (slug-ft 2 ) 

5,500 

2,210 


Example comparisons between the F-4E test results and the nonlinear 
6 DOF F-4J math model are shown in Figures 3-5. Figure 3 shows the F-4E 
at a nearly steady AOA averaging about 23 deg. The rudder trace indi¬ 
cates that the yaw damper is on. The wing rock starts as a divergent 
dutch roll but then appears to limit at approximately 40-50 deg roll 
amplitude with a $/0 ratio between 2 and 3. Also the divergence grows 
linearly, instead of exponentially as it would at constant C. During 
the large, constant-amplitude <t> and 0 oscillation, the AOA trace shows a 
small oscillation at twice the frequency of the lateral oscillation. 

Figure 4 is the F-4J simulation trimmed to a Q ■ 23 deg and disturbed 
by a lateral etick pulse. Here the dutch roll also diverges almost 
linearly for about three cycles and then becomes limited in amplitude. 
The <J>/8 ratio is about 3 during the limit cycle. An oscillation devel¬ 
ops in the AOA trace, again with a frequency double that of the lateral 
oscillation. The characteristics are remarkably similar except for the 
period of the oscillations. The higher frequency of the F-4E oscilla¬ 
tions is most likely due to a higher dynamic pressure. In Figure 3 the 
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Figure 5. Departure During Left Turn 


aircraft was pulled up to a * 23 deg within about 10 sec and then held 
at this AOA for about 15 aec. The initial Mach number waa not indi¬ 
cated. In Figure 4 the model was trimmed at an initial speed corre¬ 
sponding to o Q - 23 deg. If all other factoro were exactly the same, 
the difference In frequency shown here would be accounted for with only 
a 30 percent difference in speed. 

The computer-generated Images of Figure 5 show a comparison of the 

F-4E aircraft and the F-4J model departures from a wind-up turn to the 

left. The departure onset is quite similar. Both depart with a nose 
slice away from the turn accompanied by a roll oscillation. The F-4E 
ends up with more of a rolling departure, while the F-4J ends with more 
of a yaw departure (nose slice). This difference in final motion is 
probably influenced greatly by the differences in inertias. The F-4K 
has 24 percent higher yaw inertia and 7.3 percent higher ratio of yaw to 

roll inertia. Thus it could be expected to be more resistant to yaw 

departure and prone to roll departure. 
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Since our main interest is in modeling and establishing causal rela¬ 
tionships leading to high AOA maneuver-limiting factors rather than 
modeling post-stall gyration and spin, the foregoing match of the non¬ 
linear wing rock behavior and nose slice onset is considered to be one 
validation of the aerodynamic model. Appendix I of Part III and por¬ 
tions of the analysis and simulation to be presented later offer addi¬ 
tional validation. 

One key aspect is that the aerodynamic model demonstrates both limit 
cycle wing rock and nose slice departure without Introducing artificial 
hysteresis in aerodynamic coefficients (Reference 14), drastic changes 
in roll damping (Reference 10), or contrived elevator control asymmetry 
(Reference 15). 

3. Dynamic Analysis 

In high AOA maneuvering flight it is not uncommon for the aircraft 
to be in asymmetric (6 j 0) flight — either intentionally or uninten¬ 
tionally. Previous studies (e.g., References 5, 16-18) have shown that 
static aerodynamic cross-coupling due to sideslip can have significant 
influence on aircraft stability and dynamic characteristics and possibly 
the departure characteristics. For example, Figures 6-8 reflect the 
pitch, roll, and yaw aerodynamic momenta, respectively, for our F-4J 
model as a function of a and 8. For simplicity the aerodynamic moments 
were assumed symmetric with sideslip. Over the region 15 < a < 25 deg 
pitching moment is significantly influenced by sideslip, and rolling and 
yawing momenta due to sideslip vary greatly with AOA. This is the same 
AOA region in which the various maneuver-limiting characteristics of the 
F-4 are exhibited. Thus, in the analysis to follow, attention will be 
devoted to examining the possible influence of static aerodynamic cross¬ 
coupling on high AOA maneuver-limiting factors. Aircraft open-loop 
stability characteristics will be analyzed first, and then cloBed-loop 
piloted control aspects will be determined. 
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Figure 8. F-4J Yaw Coefficient Variation 
with a, 8 


«• Open-Loop Dynamic Characteristics 

Figure 9 presents a locus of longitudinal short-period an 1 lateral 
dutch roll root migrations with AOA in wings-level, symmetric, 1-g 
flight at an altitude of 15,000 ft. These dynamic parameters represent 
symmetric, linearized, fixed-operating-point conditions. With 
increasing AOA the dutch roll mode is seen to progress from a lightly 
damped stable mode to oscillatory divergence and, finally, to a pair of 
aperiodic divergences while the longitudinal short period remains 
relatively unchanged. 

Also shown in the figure is the AOA at which the dynamic stability 
parameter Cng^ n becomes zero. At lower AOA the parameter is positive, 
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Figure 9. F-4J Lateral-Longitudinal Short-Period Root 
Migration with otj - 0 


Indicating directional dynamic stability. At high AOA the parameter is 
negative, indicating directional dynamic instability. The AOA regions 
below and above Cngjy n ■ 0 are identified as wing-rock and nose-slice 
regions, which correlate quite well with the AOA ranges Identified with 
these characteristics in the F-4E flight test, Figure 10 (from Refer¬ 
ence 19). 

The influence of non-zero sideslip at trim on the dutch roll and 
short-period frequencies and damping 1 b reflected in Figure 11. The 
6 DOF aircraft was trimmed to l.S and 5.5 deg sideslip and linearized 
transfer function parameters obtained. These Indicate that the dutch 
roll is destabilized and the short-period stabilized by steady sideslip. 
Thus there is an apparent Interchange of damping between the two moden 
while the overall system damping remains essentially constant. 


% 

f* 


Unfortunately, comparison of dynamic characteristics predicted by 
linear frozen-point analysis in Figures 9 and 11 with the actual flight 

17 






Symbol 

Gross Weight 
(lb) 

Aim Norma! 
Load Factor 

(g) 

Symbol 

Gross Weight 
(lb) 

Aim Normal 
Load Factor 

(a) 

O 

40,500 

1,0 

★ 

35,200 

2.0 

□ 

40,000 

1.0 

* 

34,800 

3.0 

A 

39,400 

1,0 

□r 

40,000 

1,0 

T 

39,200 

1.0 

A 

39,700 

2.0 

♦ 

39,000 

1.5 

V 

30,700 

2.0 

1 

39,200 

2.0 

<? 

43,500 

1.0 
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37,100 

1.0 

A 

43,200 

1.0 

A 

36,500 

1.0 


42,700 

2.0 
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36,200 

1.0 
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42,200 
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Figure 10. F-4E Stall Approach Characteristic! 
(From Ref. 19) 
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Figure 11* F-4J 6 DOF Linearized Equation; Lateral- 
Longitudinal Root Migration with a and B 


traces of Figure 3 or nonlinear 6 DOF model traces of Figure 4 shows 
several significant differences. The wing rock linear analysis predicts 
a dutch roll divergence of increasing amplitude as a, 8, or both are 
Increased, and a longitudinal short-period with a frequency quite close 
to the dutch roll frequency, but heavily damped. The traceB of 
Figures 3 and 4 show the dutch roll divergence bounded In amplitude and 
a low-damped longitudinal oscillation at about twice the dutch roll 
frequency. 

Figure 12 shows a similar set of nonlinear time response character¬ 
istics for the open-loop 6 DOF nonlinear model trimmed to a Q ■ 21 deg 
and 0 O - 0 deg and excited with an aileron pulse. From Figure 12 it is 
apparent that the longitudinal oscillation results from rectifying the 
lateral oscillation and the sharp peaks of pitching acceleration 
coincide with sideslip passing through zero. Thus, the nonlinear-type 
behavior appears to be caused by lateral-longitudinal coupling due to 
sideslip. 

A nine-by-nine matrix of aerodynamic and kinematic terms for 
coupled, non-symmetric flight obtained from the partial derivative 
expansion of the nine equations of motion is presented in Figure 13. 
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The major derivatives and kinematic terms are identified in literal form 
and, for comparison with Figure 12, have been evaluated at a Q ■ 21 deg 
and 8 0 ■ 1.5 deg. The principal cross-coupling at this small sideslip 
angle is seen to be due to Mg, «C £ and N£. The Influences of these 

derivatives on actual vehicle dynamic response can best be seen through 

the use of eigenvectors and vector polygons. 

Figure 14 contains the eigenvectors for the dutch roll and short- 
period modes. Those show both modes to be dominated by roll, side¬ 
slip, and pitch, with the difference between the modes being primarily 

motion phasing and damping. For example, in the dutch roll, q and 8 are 
nearly 180 deg out of phase, while in the short-period they are nearly 
in phase. However, the short-period is so well damped (C a p * 0.6) that 
it would influence only the first second or two of the time responses 
shown in Figure 12. Therefore, the observed time responses must be due 
entirely to the dutch roll. 

Vector polygons for the four principal motion equations are pre¬ 
sented in Figure 15. Each polygon represents a force or moment equa¬ 
tion. The polygon's sides show the contributions of each motion varia¬ 

ble to the total acceleration. Vector lengths are a product of the 
appropriate stability derivatives and eigenvectors evaluated at the 
dutch roll frequency. The predicted motion is mostly rolling (p), with 
sideslip (8), and pitching (q) about the same magnitude, and with com¬ 
paratively little yawing (r). The polygons for o and u are so small 

that they can be neglected. The vector polygons provide the first 

insight to the cauaeti for the change in modal response characteris¬ 
tics. The p vector polygon shows phasing of the vector opposite 

to that for the <£pp vector. Thus the aerodynamic cross-coupling 
coefficient, ln effect opposes the normal aerodynamic damping 

term, £p, and therefore increases the rolling tendency of the aircraft. 
Since the aerodynamic cross-coupling coefficients Increase in magnitude 
as the trim sideslip increases, one can readily visualize that at higher 
sideslip angles, can actually become greater than £p. Then one 
could expect the model to exhibit a roll divergence for large a and 8 
perturbations. 
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Figure 15. F-4J Force Vector Folygans; 

«o “ 21 de 8. e 0 “ 1,5 de « 

The 3 vector polygon indicates that this vehicle motion derives pri¬ 
marily from the kinematic coupling of roll rate about the body center- 
line axis at the elevated AOA. 

In the r vector polygon the p vector Is missing because the aero¬ 
dynamic coefficient Np ia negligible at this flight condition. The 
vector N£a has a relatively large component in phase with N£r so that 
the aerodynamic cross-coupling from longitudinal into lateral augments 
the natural yaw damping and thus contributes to reducing body-axis 
yawing motion. All of thla contributes to the dominance of roll and 
sideslip in the Figure 12 wing rock. 

Finally the q vector polygon shows that the M^B vector is of oppo¬ 
site phase and is larger in magnitude than the damping vector derived 
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from the aerodynamic derivatives Mq + M„. Thus the sideslip-induced 
pitching moment overpowers the normal aerodynamic damping and forces the 
double frequency oscillation observed in the time traces. Mg changes 
sign with 0 (Figure 6). Since pitch inertia is considerably larger than 
roll Inertia for this aircraft, a transfer of energy from the unstable 
lateral mode into the stable longitudinal mode may be a contributing 
factor In bounding the lateral-directional mode divergence. The influ¬ 
ence of reducl.'■■ is shown in Figure 16. In this case, is small 
but not zero. Comparing these traceB with those of Figure 12 demon¬ 
strates the strong influence Mg has in bounding the lateral divergence. 

It Is concluded for our model that the roll/sideslip oscillation 
starts as a divergent dutch roll with small sideslip excursions. As 
sideslip ‘increases, coupling causes a nonlinear behavior which bounds 
the divergence. Thus F-4J wing rock, per se, is a nonlinear phenomenon 
dominated by static aerodynamic latoral-longltudinal croBs-coupling. 
Further, the vector polygons show that except for c£p the damping 
derivatives do not have a significant influence on high AOA dynamic 
characteristics. 

This analysis has also shown that one must exercise caution in the 
use or interpretation of open-loop transfer function parameters even 
when obtained from 6 DOF fixed operating point conditions. Such param¬ 
eters are limited to relatively small sideslip conditions such as onset 
to wing rock or departure and hence better reflect susceptibility to 
these maneuver-limiting factors rather than describing the resulting 
vehicle motion. 

b. Closed-Loop Dynamic Characteristics 

The cloBed-loop dynamic characteristics of any Bystem are strongly 
Influenced by transfer function numerator roots (zeros) in or near the 
range of the desired bandwidth of control. When controlling any motion 
variable with any controller, the numerator zeros of that transfer func¬ 
tion attract denominator (open-loop) poles. Hence these zeros alter 
closed-loop stability and modal response characteristics of the vehicle 
in proportion to the tightness (gain) of the loop closure. Thus certain 
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zeros can become high AOA maneuver-limiting factors in that they cause 
or influence closed-loop roll reversal, wing rock (lateral PIO), and 
departure (nose slice or roll). This observation applies to loops 
closed both by the pilot and by automatic feedback systems. The key 
numerators for piloted control are N t stlc » N $ r » and N S s tab* Tlie 
augmentation systems are of little concern due to limited SAS authori¬ 
ties and very low loop gains at high AOA and low dynamic pressure. 

Lateral stick control of roll . Under normal low AOA flight 
situations, the 3 DOF equations produce a transfer function 

. N ^stk 

A 

N S stk " + + 

A - (s + 1/T g )(s + 1/Tr)(s2 + 2c d u) d s + uj) 


<Ks) 


with values of ui^ and 5^ in the vicinity of those for the dutch roll 
mode, u> d and C d <a d . With this loop closure the zeros attract (modify) 
the dutch roll, causing the closed-loop root locus (with increasing 
feedback gain K) to take on a form typified by the sketch below. 


)w 



‘frcmcU- — 8 »tk 
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2 

As AOA Is Increased, generally decreases and can become negative 
In value. When this occurs, N | stlt 1* of the form 



where one root ie, positive and one la negative. With the high AOA shift 
In open-loop denominator roots, the root locus often takes on the form 
shown below. 



In this case the spiral root, -1/T gl Is driven toward the right 
half-plane (RHP) aero and a first-order Instability (divergence) re¬ 
sults. The rate of divergence depends upon how far the zero lies in the 
RHP and how tightly the loop is closed. 

Migration of the F-4J open-loop zeros with AOA at zero sideslip 
Is reflected by the solid lines In Figure 17, which indicate an almost 
2 

linear decrease In “a (at constant damping) with increase in AOA up to 
• 2 

a - 18 deg. Above this angle, la negative and the RHP root location 
is quite AOA-sensitive. Also noted on the figure are the regions asso¬ 
ciated with wing rock and nose slice from Figure 10. Wing rock occurs 
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in the region where is emallj nose slice occurs in the region where 
1/T,^ is large and negative. 

2 

Actually, negative reflects steady-state roll opposite to that 
commanded by stick deflection. This is generally due to large adverse 
aileron yaw, since (after Reference 20) 



where 6^^ “ 8<5 fl + 11.66 ep . At high AOA, Ng is small, iCg is large 
negative, and N 6 stk determines the sign of «<)>• Adverse (negative) 
then tends to produce negative and the RHP zero. In physical terms, 
adverse yaw produces sideslip and, in turn, roll reversal due to £g0. 
It thus becomes apparent that the nose slice region identified in Fig¬ 
ure 17 may be related to strong adverse aileron yaw and the magnitude of 
the RHP zero. 

The dashed lines of Figure 17 show the values of ^ and 2?^^ ( 0 r 
l/T^, l/T,^) for e o " 5,5 de 8* Wherea3 increasing AOA at 0 ■ 0 causes 
u>? to decrease steadily and finally become negative, when S M the 
^ 2 

aerodynamic cross-coupling causes to stay positive and close to the 
dutch roll mode but the damping ratio, becomes negative (see sketch 
below). 


Lateral - Directional 
Root Locations for: 



0o = 5.5° 
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The promlxity of u>^ and means low residue for the dutch roll mode 
In roll control with stick when $ 1 * 0. Thus sideslip should reduce the 
roll content of the unstable dutch roll mode; this is consistent with 
the observed bounding of the roll oscillation divergence In Figures 4 
and 16. 

Rudder control of roll . When the region of adverse aileron yaw Is 
reached, roll control is maintained via rudder. This produces yaw and 
sideslip-induced roll in the same direction and thus much more rapid 
turning than just yawing moment and side force would produce. The 
rudder-to-roll numerator, N$ r » of the form 



One root is always positive (RHP) as shown in Figure 18. The RHP zero 
results from adverse roll due to £fi r . It is apparent from Figure 18 
that the zeros of are not greatly influenced by either AOA or side¬ 
slip. 

Again this control structure portends closed-loop Instability since 
some pole must be driven toward the RHP zero. However, these zeros lie 
at high frequencies relative to the low-frequency bandwidth normally 
achieved with rudder control via the pilot's legs. Rudder-to-roll con¬ 
trol handling qualities are normally dominated by the low-frequency 
spiral and roll subsidence modes, and at high AOA these remain stable 
and well behaved. While this type of control is "unnatural," it does 
not result in roll stability problems at low gain closures. 
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Figure 18. Nj r Root Migration with A0A{ F-4J 

Stabilator control of pitch . Cloaed-loop control of pitch attitude 
becomes pertinent to maneuver-limiting if and when the numerator for 
stabilator control of pitch, N fi sta ^» zeros in the RHP. At zero 

sideslip the transfer function is of the form 

0(s) _ N6 atab 

-t— 



A - [s 2 + 2CpW p s + ui$][s 2 + 2c Bp aJ ap s + u» 2 p ] 
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On the F-4J, 1 /Tq 2 decreases with increasing AOA and finally couples 
with 1 /Tq^ to form a complex pair which almost cancels the longitudinal 
phugoid mode. Thus these zeros remain well behaved. However, when 
there is sideslip so that static aerodynamic cross-coupling becomes 
significant, the numerator contains additional zeros related to lateral- 
directional modes. 

It was discovered (Reference 5) that such coupling in the A-7 air¬ 
craft produced an RHP zero, l/Tg^, and subsequent pilot control of pitch 
attitude via the stabilator caused a closed-loop directional divergence 
which has the appearance of a nose slice. A similar situation exists 
in the F-4J aircraft. Figure 19 presents pole and zero locations for 
a Q - 23 deg and 6 0 - 5.3 deg. The coupled 6 DOF equations (Figure 19a) 
produce two real zeros (one RHP, one LHP), one complex pair which nearly 
cancels the phugoid, and one higher-frequency complex pair. 

The source of the low-frequency zeros can be determined by elimin¬ 
ating the velocity equation, which should reduce the phugoid mode to 
a first-order pole at the origin and eliminate the speed numerator zero, 
l/Tg^. It may be observed in Figure 19b that the only change is in the 
phugoid dipole pair which became first-order as expected. 

Next, as in Reference 5, the cross-coupling terms and are set 
to zero in the 6 DOF equations to Identify the lateral-directional 
modes. This is Bhown in Figure 19c where the pole-zero cancellations 
show the two real zeroB to be related to the roll subsidence, 1/T R , and 
spiral, 1/T # , modes and the higher-frequeucy complex zeros to be related 
to the dutch roll. 

Figure 20 is a system survey for closure of the pitch attitude loop 
for the augmented (SAS on) airframe at this same trim condition. The 
lower figure is a root Iocub and the upper figure is a Bode-siggie 
(Reference 20). The root locus shows the migration of roots from the 


*The slight difference in airframe dynamic parameters between 
Figures 20 and 11 indicates the very minor influence of the SASs at high 
AOA and low dynamic pressure conditions. 
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open-loop poles toward the open-loop zeros as gain Is Increased. It 
provides Information concerning closed-loop frequency and damping as a 
function of gain but does not provide Information concerning sensitivity 
to gain change. The Bode-slggle presents the open-loop transfer func¬ 
tion amplitude and phase as functions of frequency where the amplitude 
is evaluated for s ■ Jw (Bode) and s ■ o (Biggie). The -f. symbols iden¬ 
tify migration of complex LHP closed-loop pole frequency with increasing 
gain (Bode-root locus, Reference 20) and the 0 symbols similar migra¬ 
tion of closed-loop pries in the RHP. Closed-loop first-order roots 
lie at the intersection of the gain closure lines and the a amplitude 
curves. Thus these survey plots allow Identification of closed-loop 
roots at a specific gain ”ilue, and sensitivity of root value to changes 
in gain. 

In order to provide any pitch attitude control whatsoever, the loop 
gain must he sufficiently high that the closure line crosses the Bode 
plot below the -20 dB low-frequency asymptote, i.e., with a gain 20 dB 
or greater. Generally the crossover will be at or above the short- 

period frequency - again in the vicinity of -30 to -40 dB. The root 

locus plot has tick marks showing closed-loop root locations for gains 
of 20 and 32 dB. This shows that the root moving toward 1/Tg 3 has 
already reached the RHP, and divergence, with the gain at 20 dB. By 
32 dB the firat-order divergence has a time constant of 4 sec. Thus any 
attempted control of pitch attitude reoults in a divergence. This 
divergence is due to lateral-longitudinal coupling and is evidenced in a 
lateral motion. 

Figure 21 shows values for the real zeros 1 and 1 /Tr at 0 O ■ 
5.5 deg and <* 0 from 21 to 25 deg. These indicate that at non-zero 
sideslip and A0A Increasing above 20 deg 1 /Tq 3 rapidly moves to negative 
values which portend closed-loop first-order divergence time constants 
as low as 1 sec. Again this coincides with the nose slice region 
observed in flight test. 
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The effective transfer function for each control loop becomes 


. "S.t.b-vC.X.k’ 
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where N 6 -ta b°stk i8 the numerator term due to cross-coupling effects. 
Both the 0 and $ numerators are modified by the same coupling numerator; 
the extent of the influence depends upon how tightly the other loop is 
being closed. This Influence can be determined by treating the 
effective numerator as a closed-loop system and identifying the 
migration of numerator roots as the loop gain (Y^ or Yg) is varied, 
e.g., the roots of the 6 numerator are obtained from 


Nfi stab " Nfi stab + Y *%tab«stk 

, ,. yS;t.bt.tb 


The uncoupled numerator roots become the poles of this system and 
the coupling numerator roots are the zeros. The root migration is 
controlled by the pilot gain, Y^. For the Figure 20 example case of 
a “ 23 deg, B 0 " 5.5 deg, the transfer function is 


*• 0029 (**1.35X1.69)[ .993; .102 m 


The reader is referred to page 163 et seq. of Reference 20 for a 
thorough development and treatise of multiloop control systems analysis. 




where first-order factors (s + a) are represented by (a) and second- 
order factors [s 2 + 2cu n s + usjj] by [C5 u n ] • Assuming Y^ to be a simple 
positive gain, this produces a root locus of the form shown in the 
sketch below. Thus pilot control of the roll attitude loop via aileron 
and spoiler causes the RHP zero of the pitch numerator to move further 
into the RHP because the root of the coupling numerator lies beyond 
I/T03. 
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Similarly, pilot control of pitch attitude via etabllator alters the 
roll numerator as 

^stab^stk “.0029(-1*35)<1.69)[.993; .102]Y e 

—" Turre t v^r , "' ny T F T ggBTT:w - n . w. t ray 

4 stk 

which, for Yg assumed a simple negative gain, produces a root locus of 
the form sketched below. This results in the RHP complex pair being 
driven to the real axis with one root eventually moving further into the 
RHP while the other moves into the LHP. 













. 0041 . 7 ( 1 . 80 )(- 3 . 04 ) [ . 180 ; . 148 ][. 277 ; 1 . 20 ] 


Use of rudder to roll modifies the pitch numerator In much the same 
manner «• doeB use of aileron and spoiler. This can be seen by 
comparing values of the coupling numerator terms. On the other hand, 
pitch control has little Influence on the rudder roll numerator because 
nearly all numerator and denominator roots cancel. Thus It appears that 
this control strategy produces the most consistent closed-loop lateral- 
directional dynamic properties at high AOA and non-zero sideslip, which 
also might be a reason pilots do not mind using this technique. 


9ummsry . This analysis has shown several areas in which control- 
numerator dynamic parameters can become maneuver-limiting factors for 
the F-4J. These limits are generally associated with zeros migrating to 
the RHP at high AOA In asymmetric (0 + 0) flight. 

The RHP zero of N l Btk 1* due to adverse "aileron" yaw (Nj^ n®g«“ 
tlve) and rasults in roll reversal. The attendant sideslip causes the 
♦ zeros to shift toward the dutch roll poles. This proximity reduces 
the roll contribution in the dutch roll mode and produces a nonlinear 
bounding of dutch roll (or wing rock) oscillation amplitude. 

One RHP zero exists In whether the aircraft is in asymmetric 
flight or not. However, that root is at a frequency sufficiently above 
the region in which rudder control is generally exerted that it does not 
pose a closed-loop stability problem. 
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An RHP zero can occur in N l ata ^ due to aerodynamic cross-coupling 
with non-zero sideslip. This zero results in a first-order divergence 
mode if stabilator control of pitch attitude is attempted. 

Asymmetric flight also produces lateral-longitudinal control cou¬ 
pling that further aggravates these undesirable RHP zero locations of 
the key 6 and numerators. 

Thus the aerodynamic cross-coupling associated with sideslip tends 
to dominate the F-4J high AOA flying qualities. Values of the cross- 
coupling derivatives evaluated at B 0 “ 5.5 deg and a range of AOA are 
shown in Figure 22. Comparing these with the root migrations in 
Figures 9 and 11 demonstrates their significance. The regions where 
large shifts in dynamic parameters occur coincide with large magnitude 
of cross-coupling derivatives. Also, it should be noted that the 
magnitude of the cross-coupling derivatives is directly proportional to 
the magnitude of sideslip. 

The foregoing analytic results must be viewed with some caution 
since it was shown that linearized, fixed-operatlng-polnt dynamic 
parameters do not necessarily reflect the actual open-loop dynamics of 
the nonlinear airframe. On the other hand, the apparent correlation 
batween bounaing of dutch roll divergence end movement of zeros with 



Figure 22. F-4J Static Aerodynamic Cross-Coupling 
Derivatives; 0 O « 5.5 deg 
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sideslip, and Che AOA region for nose slice observed in actual flight 
with the region where the RHP zeros of N I s tab an< * K ^stk are * ar * e ls 
quite encouraging. 

B. F-14A ANALYSIS 

The F-14A configuration does not have wing leading-edge slats. This 
aircraft has variable wing sweep; at M < 0.5 the wings are in the full 
forward (A ■ 22 deg) position. The unaugmented airframe ls classified 
as departure resistant, but it can be departed (References 21 and 22). 
It remains laterally stable and longitudinally controllable at full aft 
stick (AOA approximately 50 deg with wings forward). It haB a mild 
dutch roll instability starting at about 15 deg AOA. The oscillations 
are aggravated if the pilot attempts to oppose the motion with either 
lateral stick or rudder. The dutch roll mode becomes stable again at 
about 23-25 deg AOA. Roll reversal due to adverse yaw from the aileron 
and differential horizontal stabilizer also starts at about 18 deg 
AOA. Prolonged deflection of the lateral stick at a > 2u deg will 
result in a departure, characterized by a rapid increase in “adverse" 
yaw rate, which may develop into a series of uncommandud rolls with the 
appearance of snap rolls. At low speed, departure can occur at AOA as 
low as 21 deg. The departure characteristics are a function of wing 
sweep, c.g. location, and Mach number. The airplane will spin and has a 
non-racoverable flat spin mode. 

1. Analytic Model 

The F-14A aerodynamic model was derived from several sources (Refer¬ 
ences 23-27). The data base was taken from the NASA/Langlev Differen¬ 
tial Maneuvering Simulator (DMS) report of Reference 23. This fixed- 
base piloted simulation was considered by Grumman test pilots to provide 
a valid representation of high AOA dynamic characteristics observed in 
early flight tests. However, time traces obtained from ther<o same aero¬ 
dynamic data did not match available flight test records. Therefore 
the date were modified on the basis of Information obtrined from other 








references! discussions with F-14 aerodynamlcists at Grumman, and 
finally, as necessary to match available flight test data. 

The resulting aerodynamic model, together with supportive documenta¬ 
tion for each aerodynamic coefficient, Is detailed in Part III, Appendix 
II. The data caver the range! 

0 < a < 55 deg; -20 < 6 < +20 deg 


at low Mach (0.2). At this low speed the variable-sweep wing is 
nominally in its full forward (A - 22 deg) position, ho wing sweep 
Influence is not included. 

Based upon the F-4J analysis, the key high AOA aerodynamic moments 
are Ck/OjB), ^(0,0) and C m (a,g). The roll moment (Figure 2.3) is seen 
to remain essentially constant for a > 15 deg and 0 < 10 deg, as opposed 
to the F-4J's significant decrease in G^(a,6) for a > 15 deg (Figure 7). 
The yaw moment (Figure 24) is Bimilar to the F-4J's (Figure 8) for 
a < 17 deg, but for a > 17 deg both the slopes with respect to a and the 
maximum values with respect to 0 era approximately half those of the 
F-4J. The combination of less negative C n ^ and larger negative Cig pro¬ 
duces positive values of the dynamic stability parameter Gngd yn through¬ 
out the 55 deg AOA range. 

The pitch moment with sideslip (Figure 25) is also similar to that 
of the P-4J ovor the range 10 < n < 25 deg, where it is negative. How¬ 
ever, above a ■ 35 deg, G m ^ becomes positive. 

Weight, c.g., and moments and product of inartia for the clean (no 
external stores) aircraft were provided by Grumman. These are shown in 
Table 2 with F-4J values for comparison. All Inertias are in the 
fuselage centerline reference axis system. Aircraft dimensions are also 
presented for comparison. 


43 




















TABLE 2 


COMPARISON OF F-14A AND F-4J WEIGHTS, 
INERTIAS, AND DIMENSIONS 


Parameter 

F-14A 

F-4J 

W (lb) 

46,950 

37,000 

c.g. (X c) 

16.0 

29.3 

I x (alug-ft 2 ) 

58,950 

23,850 

I y (slug-ft 2 ) 

225,600 

127,400 

I z (slug-ft 2 ) 

285,000 

146,000 

I xz (slug-ft 2 ) 

-3030 

2210 

A (deg) 

22 

45 

s (ft 2 ) 

565 

530 

b (ft) 

64.1 

38.67 

c (ft) 

9.8 

16.04 


The inertia differences between the vehicles are large. Although 
the F-14A ia some 27 percent heavier than the F-4J with comparable 
fuel loading, the inertias range from 77 to 147 percent higher and the 
product of inertia is of opposite sign. This has considerable influence 
on the ratios of aerodynamic to kinematic terms in the model and greatly 
influences coupling effects, as will he shown later. 

Aircraft control is exerted through a horizontal stabilizer (6 ata b), 
twin rudders (6 r ), spoilers and differential horizontal (<$ Q ) sur¬ 
face deflections. Simplified block diagrams of the flight control and 
augmentation systems are shown in Figures 26-28. The longitudinal sys¬ 
tem contains a feel system (spring cartridge, variable gearing ratio, 
bobweights), a limited-authority pitch rate SAS, and a surface actuator. 
However, the analysis and simulation are based on a simple fixed stlck- 
to-Burfaoe gain and pitch rate SAS. 
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Figure 28. F-14A Directional Control System 


The lateral system (Figure 27) has a roll rate command augmentation 
system (CAS) in parallel with a mechanical stick-surface link. At high 
AOA the spoilers are Ineffective and all roll control is obtained 
through differential horizontal, The roll rate CAS operates only 

through the differential horizontal, and the electrical path gains are 
varied with AOA such that the CAS command and feedback fades to zero 
when 20 < a < 31 deg and, simultaneously, the electric stick-Burface 
path fades in to cancel all but ±2 deg of mechanical path surface 
command. Thus, for a > 31 deg there Is no roll rate damping augmenta¬ 
tion and the pilot can only command ±2 deg of differential horizontal 
surface. 

The rudder system (Figure 28) contains a fixed gain mechanical link 
from pedals to surface actuators, a conventional yaw damper with yaw 
rate and lateral acceleration feedback, and a lateral stlck-to-rudder 
crossfeed (SRI). At a > 10 deg the yaw SAS is cut out and the SRI faded 















in. The SRI signal, passed through a first-order lag filter (break 
frequency 8 rad/sec), reaches full gain at a > 20 deg. Above a ■ 20 deg 
full lateral stick deflection produces 19 deg of rudder deflection. 

Thus for the AOA region of interest, a > 20 deg, the airframe is 
essentially unaugmented in the lnteral-directional axes and the control 
stick effectively moves the stabilator for pitching and twin rudders for 
rolling. 

2. Model Validation 

F-14A flight test data covering the appropriate AOA region were 
provided by the aircraft manufacturer (Reference 28). Unfortunately, 
the data were extracted from tests of three aircraft having slightly 
different weights, inertias, flight control systems, and flight test 
instrumentation. Since the flight test goals Involved flight control 
and leading edge flap system development rather chan identification of 
basic airframe dynamics, most runs were contaminated by SAS effects, 
pilot control inputs, or configuration influences. The model validation 
therefore was based on general response matching of a few short time 
"windows" of bare-airframc response or rel*i-*aly uncontaminated traces, 
plus observation that the sequence of dynamic characteristic changes 
with increasing AOA described previously was indeed achieved. 

Figure 29 presents results of a static longitudinal validation check 
based upon a 1-g stall in which AOA is increased steadily from 10 deg to 
almost 40 deg with a ramp horizontal stabilizer input. The circled 
points reflect the excellent agreement for horizontal stabilizer inputs 
required to achieve the noted AOA trim condition at 35,000 ft altitude 
with the aerodynamic model. 

Figure 30 is a short time-slice of the divergent dutch rolx at 
approximately 15 deg AOA. Motion 1 b excited by a lateral stick doublet 
which, through the SRI, produces both differential horizontal and rudder 
doublets. Both the yaw and roll rate SA3s are then turned off, and the 
subsequent oscillation represents the "bare" airframe in the lateral- 
directional sense. As with the F-4J, the envelope of the oscillation 


49 







35000 

33000 

[31000 

29000 


p j. rr r \r, j-jj a j-v r a r^'j JT-rf . -M B 

PiliSIOlEiTtE 

? [feUr'i-ffcl'H: U--4ij uph. 

rH-rt ! .**}•:■ I tf t* 1 i-t-t-h-Frt-ht*r 


a3Xoasi:r4:pqj.t.EXtjJ: 


.kliitiil 





Figure 29. F-14A Longitudinal Static Trim Comparison 
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seems to grow linearly rather than exponentially, which indicates 
nonlinear damping. AIbo, the AOA trace shows an oscillation at twice 
the dutch roll frequency. The motion on the 6 sta j, trace indicates the 
pitch SAS was on and responding to the longitudinal oscillation. 

Figure 31 presents time traces for the F-14 aerodynamic model 
trimmed at 15 deg AOA and excited by aileron and rudder doublets similar 
to Figure 30. This also shows an almost linear, rather than 
exponential, growth in dutch roll amplitude. The frequency and cycles 
to double amplitude are the same as in Figure 30. The small-amplitude 
AOA oscillation also is precisely twice the frequency of the dutch roll. 

Figure 32 presents another short time-slice in which AOA is 

maintained at 18 ± 2 deg. In this angle of attack region the SAS gain 

schedule shows that the yaw SAS should be off and the roll SAS on. 

However, comparison of the yaw rate and rudder position traces indicates 
that the rudder is correlated with yaw rate at approximately the 
magnitude and phase angle that would result from the yaw damper. Since 
the rudder pedal and lateral stick are both constant end very small, the 
aircraft response is essentially that of the controls-flxed augmented 
airframe. This dutch roll oscillation is essentially the same, 
indicating that the yaw SAS is quite ineffective at this AOA and 

airspeed. 

Other comparisons which can be made among Figures 30, 31, and 32 are 
summarized below: 


Parameter 

F- 

Figure 32 

■14A 

Figure 30 

Period of oscillation (Bee) 

4.5 

4.5 

Cycles to double 

2 

2 

$/3 (deg/deg) 

2 

* 

p/r 

25 

** 

a (deg) 

18 ± 2 

15 


Model 
Figure 31 

4.5 
2 
2 

22.5 
15 


@ trace not available. 


r trace contaminated. 
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Taken together, the foregoing Is considered an acceptable match between 
actual and simulated airframe dynamics in the region of 15 deg AOA where 
the dutch roll divergence and wing rock are most pronounced. 

Figure 33 shows a pullup initiated from wings-level flight in which 
AOA is steadily increased from 13 to 38 deg at about 0.75 deg/sec. A 
lateral stick oscillation is introduced at about 15 deg AOA. The stick 
is then centered and the SASs turned off at about 18 deg AOA. A small 
rudder ramp is then introduced, possibly to counter the Blow roll-off in 
bank angle. Simulation traces for a matching pullup with similar stick 
and rudder inputs are shown in Figure 34. Both Figures 33 and 34 show 
a divergent dutch roll. However, the amplitude appears to increase 
linearly rather than exponentially, again reflecting nonlinear damping. 
Both sets of traces show the period of oscillation to Increase from 
about 4.5 sec to over 5.5 sec as AOA approaches 30 deg. Again there is 
good agreement between the actual flight and simulation responses. 

3. Dynamic Analysis 

The principal shortcomings of the F-14A bare-airframe flight char¬ 
acteristics at high AOA are wing rock (unstable dutch roll) and roll 
reversal. Otherwise the airplane is quite docile. It is departure- 
resistant, but can be departed if lateral control deflection iB sus¬ 
tained without large accompanying rudder. Thus, the following discus¬ 
sion will be devoted to identifying potential factors underlying these 
key differences between the F-14A and F-4J. 

a. Open-Loop Dynamic Gharactaristlcn 

A Burvey plot of root migration with angle of attack for the range 
15 < oi < 30 deg and 8 * 0 is shown in Figure 35. This indicates excel¬ 
lent lateral-directional agreement with the previously reported high AOA 
flight characteristics and with the flight traces of Figures 30, 31, and 
33. The dutch roll mode becomes divergent at just under 15 deg AOA and 
stays unstable until approximately 22-23 deg AOA. This is the region of 
reported wing rock. Above 23 deg AOA the dutch roll is stable and well 
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Figure 34. F-14A Model Simulation of 

Flight 236-1 Pullup 
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Figure 35. F-14A Root Migrations with AOA 


damped but the coupled roll-subsidence and spiral becomes oscillatory. 
The AOA at which ■ 0 is not shown on this plot because it remains 
positive (as shown in Figure 36). 

The longitudinal short-period mode exhibits relatively low frequency 
and high damping at 15 deg AOA. Frequency then increases and damping 
decreases with increasing AOA until the dutch roll and short-period are 
about equal at 30 deg AOA. Note that this is not in agreement with the 
flight traces for 15 deg AOA (Figure 30), which showed a short-period 
frequency twice that of the dutch roll. 

Referring back to Figure 25, 15 deg AOA is where Cmg is a maximum. 
Thus, cross-coupling through Mg may again create pitch frequency cou¬ 
pling and contribute to the nonlinear dutch roll damping, as was the 
case with the F-4J. However for 6 + 0 the shift in the Figure 35 root 
locations is so small as to be negligible, and this is not the same as 
with the F-4J. 


58 






1.8 

1.6 
1.4 

N /3dyn r 

1.0 
,8 
.6 
.4 
.2 
0 

Figure 36. N Bdyft Variation with AOAi F-14A 

Eigenvectors for the coupled dutch roll and short-period modes 
of the 6 DOF airframe at <x Q - 20 deg and 0 O - 1.5 dag are shown in 
Figure 37. There is essentially no longitudinal coupling into the dutch 
roll mode and only alight lateral coupling (mainly rolling motion) into 
the short-period mode. These eigenvectors are noticeably different from 
those for the F-4J (see Figure 14) at a similar flight condition. 

Vector polygons for the dutch roll mode at this flight condition 
are shown in Figure 38. The p equation indicates that X B and <C ' com¬ 
pletely dominate rolling motion. All damping derives from X'; there 
is no cross-coupling influence from All terms in the * equation 

have been doubled to make this polygon legible. The relative size 
of the p end r polygons indicates that in this AOA region static and 
dynamic yawing force end moment terms are negligible in the dutch roll 
mode. However, the yawing moment that doea exist derives primarily from 
Np. The 4 equation has significant magnitude compared to the P equa¬ 
tion. The key is the dominance of the Mo 101 term which, as with the 
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Figure 38. F-14A Vector Polygons for Dutch Roll Mode; 

a 0 - 20 deg, 3 0 “ 1.5 deg 







F-4J, indicates that coupling from the lateral into the longitudinal 
produces the frequency doubling effect. 

Based upon the above argument, it 1 b concluded that the F-14A wing 
rock 1 b due to a divergent dutch roll which derives from very low 
damping in roll, jCp, rather than lateral-longitudinal cross-coupling 
(since the open-loop lateral-directional dynamic characteristics are 
relatively insensitive to sideslip). The principal static aerodynamic 
cross-coupling is a transfer of energy from rolling-sideslippiug motion 
into pitch motion and this, as with the F-4J, results in frequency 
doubling into pitch and may contribute to amplitude bounding of the wing 
rock. 

b. Closed-Loop Dynamic Characteristics 

Lateral stick control of roll . Figure 39 presents the key 
parameters for the numerator The lower plot indicates that the 

movement of one zero into the RHP with increasing AOA is not as rapid 
for this airframe as for the F-4J (see Figure 17). However, departure 
can occur where one of the real roots is large negative. This plot also 
demonstrates that the roll numerator zeros are quite insensitive to 
sideslip when compared to the F-4J. 

Rudder control of roll . Same characteristics as for the F-4J. 

Stabilator control of pitch . The potential for divergence lies 
with the RHP zero, I/Tq^, in Figure 40 is a plot of l/Tg^ at 

S 0 “ 4 deg over a range of AOA. Note that the parameter scale factor is 
an order of magnitude less than the one used for the F-4J in Figure 21, 
and the largest value (approximately -0.01) is two orders of magnitude 
less than for the F-4J (equal to -1.0). Thus, pilot control of pitch 
attitude should have no noticeable coupling into the lateral-directional 
modes for the F-14A. 
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Figure 40. N 6 ata t First-Order Root Migration with AQA, F-14A 


Coupling numerator . The lead coefficient (gain term) for the F-14A 
coupling numerator, N fi atat) t D * is, for all practical purpoaea, zero. 
Thua closure of either loop (pitch or roll) has essentially no influence 
on the other. 

Summary . The F-14A exhibits somewhat better open-loop static and 
dynamic stability than does the F-4J. It also has considerably less 
lateral-longitudinal static cross-coupling in asymmetric flight, 
and hence little change or degradation in dynamic parameters with 
sideslip. The static aerodynamic cross-coupling derivative values at 
S Q ■ 4 deg and 15 < o < 30 deg are presented in Figure 41. Comparison 



Figure 41. F-14A Static Aerodynamic Cross-Coupling 

Derivatives 
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of these values with those for the F-4J (Figure 22) shows a drastic 
difference between the two aircraft. 

The most significant maneuver-limiting factor for the F-14A appears 
to be roll reversal associated with movement of one root of the roll 
numerator, N$ D> into the RHP as AOA is increased. The F-14A exhibits a 
more gradual shift of this root with AOA than does the F-4J. 
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section hi 


AEEODTNAHIC CROSS-COUPLING INFLUENCES 


The previous section showed a significant migration of denominator 
and numerator roots for the F-4J due to aerodynamic cross-coupling at 
0 O ■ 5.5 deg. In this section literal analytic expressions are devel¬ 
oped to examine in greater detail the influences of specific cross¬ 
coupling coefficients on the airframe characteristic (open-loop) modes 
and the key piloted control numerators for roll due to lateral stick and 
pitch due to stabllator deflections. 

With lateral-longitudinal cross-coupling, at least 5 DOF equations 
of motion (e.g., Figure 13) must be used for analysis. As a result, the 
wealth of existing knowledge on the dynamics of vehicles in three (or 
fewer) degrees of freedom is not adequate, and qualitative insight into 
the motions of the aircraft is sacrificed. 

The Important cross-coupling terms in the equations of motion iden¬ 
tified empirically by Porter and Loomis (Reference 29) and Johnston and 
Hogge (Reference 5) are both aerodynamic and kinematic. Effects of 
cross-coupling on the characteristic aircraft modes were evaluated by 
Hamel (Reference 16) by applying servo analysis techniques to a 5 DOF 
aircraft mathematical model. 

The approach applied next will follow and expand upon that utilized 
by Hamel. 

A. CHARACTERISTICS MODES 

Analysis of an uncoupled (0 = 0) airframe is relatively straight¬ 
forward. Approximations to the classical 3 DOF lateral or longitudinal 
expressions can be developed for a wide range of airframe and flight 
conditions (e.g., Reference 20). 

However, sideslip characteristically introduces cross-coupling, both 
kinematic and aerodynamic. Then the only clear simplification one can 
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apply Is the assumption of constant forward velocity, thereby removing 
the velocity equation and associated terms (see Figure 13). This can 
be justified since forward velocity changes comparatively slowly. No 
significant coupling occurs due to velocity, and its effect is primarily 
on the low-frequency characteristics (i.e., the phugoid mode). The 
three kinematic equations of Figure 13 can be resolved and substituted 
in the Z, M, and Y equations to reduce the number of equations to five. 

Hamel takes five-degree-of-freedom equations of motion "using body 
axes which have the x-axis initially aligned with the flight path"; 
from the context we interpret that for initial sideslip, the x-axis is 
aligned with the projection of the flight path vector onto the plane of 
symmetry. In expanding the characteristic determinant he separates out 
the effects of the aerodynamic coupling terms JCq, N„ and Mg (neglect¬ 
ing any other coupling between the longitudinal set and the lateral- 
directional set) thus: 

A " A long^lat + K (0 o ) A coupl 

The uncoupled (3 Q S 0) denominator terms are 

Along * a 2 + 2c sp ui ap s + u>l p 

Alat " (• + !/*.)(■ + l/T R )(s2 + 2; d u> d s + s»J) 

Hamel defines the classical longitudinal (short period) and lateral- 
directional (spiral, roll subsidence, dutch roll) by approximate factors 
(Reference 20). The coupling component, K(0 o )A C oupl' can ** manipulated 
to represent a sum of two terms: 

K(0o>A coupl - ^B 0 s(s + 1/T 0 )(.2 + 2c 0 u> 0 s + ttg) 

+ M 0 N'B o (s + 1/TiXs + 1/T 2 ) 

The approximate factors for these coupling terms are reproduced (in the 
notation of this report) in Table 3. 
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TABLE 3. COUPLING NUMERATOR FACTORS 


L/T 0 - -M q 

2Co“>o " 

<*>2 - Ng + YgNj. - ; + Yg£j) 

1/T], " “£p + (^a/ N a)(^p ~ s/^o) 

i/x 2 - t i 8 /[u 0 (A^/n'Kn; - £;)] 

The roots of the coupling component are evaluated by applying servo 
analysis techniques for obtaining the roots of a characteristic equa¬ 
tion, setting it into the form 

i + M g(V*«H s + 1 /T lX» + l/T }) m Q 

s(s + f/T 0 )[s z + _ 2C p u 0 8 + *2] 

Figure 42 is a root locus plot for this equation. In this instance 
5 0 > 1 and u> 2 < 0. The "cloped-loop” roots for this case are 

K(0 o )A C oupl " J^M 8 + *404)(s + .47) 

x [a 2 + 2(-.48)(1.87)s + (1.87) 2 j 

This Bame technique can be applied on a larger scale to find the 
roots for the coupled charactf '.atic modes, i.e., the roots of the rela¬ 
tionship 








Figure 42. Migration of Roots of K(B 0 )A C0up i with 
Increasing M & (« 0 - 21 deg, B 0 ■ 5.5 deg) 


Figure 43 is a root locus plot for this expression showing the migration 
of ia 8p and ui d with increasing X^o* Note that !•/*„ and 1 / t r are 
essentially cancelled by zeros obtained from the approximate factors and 
the preceding "closure." 

The actual 6 DOF roots of the F-4J a Q • 21 deg, (J 0 - 5.5 deg case 
are shown by the dashed lines and & symbols in Figure 43. There is 
fair agreement between the simplified and exact loci for the dutch roll 
but not for the other modes. The problem, it appears, lies with the 
location of the first-order coupling zeroB. If these were to lie closer 
to the origin, the ui sp and w d loci would be rotated counterclockwise and 
1/T a and 1/T R would be driven to the right. It is not surprising that 
the approximate factors employed here do not precisely predict the F-4J 
root shift since they were developed in Reference 16 for an airframe 
configuration considerably different from the F-4J. However, they do 
demonstrate the trends and show specifically that: 
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A Actual 6 DOF Roots 



Figure 43. Migration of Roots of Characteristic 
Modes ub Cross-Coupling Increases; 
a Q - 21 deg, 0 O - 5.5 deg 


* The location of the roots of K(B 0 )& c0upl (Fig¬ 
ure 42) and hence the direct Io n of the root shift 
in Figure 43, is controlled by the sign and mag¬ 
nitude of M p (Ng/£'). 

• The magnitude of root shift in Figure 43 is con¬ 
trolled by the sign and magnitude of cC. £ 0 O . 

B. n| NUMERATOR 

The first step in identifying the influences of cross-coupling terms 
on the key numerator for lateral piloting control, N ? st i c * require elimi¬ 
nation of all insignificant terms. This is accomplished by a series 
of simplifying steps. The 6 DOF characteristic poles and the zeros 
for the numerator N t atk «re shown in Figure 44, again for the F-4J at 
« c - 21 deg, 8 q ■ 5.5 deg. Figure 44a represents a case for which the 
cross-coupling terms are zero. At this high AOA the uncoupled numerator 
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Figure 44. N^/A Pole-Zero Locutions, Coupled vs. Uncoupled, 
F-4J, a Q - 21 deg, - 5.5 deg 


consists of two real roots, 1/T^ end l/T^, non-minimum phase and 
minimum phase, respectively; and two complex pairs, and u>$ p , which 
exactly cancel the uncoupled longitudinal poles. With croBs-coupling 
(Figure 44b) the first-order zeros couple into a complex pair located 
near the divergent dutch roll mode. In addition, the numerator counter¬ 
part to u 8 p, m{ p , has moved apart from the short-period poles. The 
effects of coupling on the phugoid numerator term, are neglible (see 
Figure 44b), This indicates that the velocity equation of Figure 13 may 
be eliminated. The phugoid then reduces to a first-order pole-zero pair 
at the origin, but all other roots are affected only slightly (Figure 
44c). However, even with the Vj. terms removed from Figure 13, and 
recognizing that the equations for the Euler angles i|i and 9 contribute 
no significant kinematic coupling and can be deleted, a rather unwieldy 
six-by-six matrix still remains. 

The matrix can be further reduced by use of time vectors evaluated 
for the particular mode of Interest. Inspection of Figure 44b shows 
that the N t SC | ( zeros of greatest concern are the complex pair which lie 
in the right half-plane, near the already divergent dutch roll poles. 
Therefore further analysis of n| ^ concentrate on this complex 
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pair, and on the effects of cross-coupling in driving them off the real 
axis, where Figure 44a places the corresponding zeros of the uncoupled 
response. 

For denominator roots the relative magnitude and phase of each 
motion variable in the mode corresponding to any particular real pole 
or pair of complex poles can be determined. The result can be expressed 
numerically as an eigenvector or graphically as phaeors. In order to 
determine such characteristics for transfer function numerator roots, a 
pole is driven into each zero tu be examined. The poles are driven Into 
the zeros by closing a tight feedback loop, as shown in Figure 45. In 
Figure 44c the real poles are driven into the right half-plane pair of 
complex zeros. Actually the (p/s)6 fltlt transfer function is utilized 
here rather than $/6 at k since the two are almost identical (compare 
Figures 44c and 45). 

Figure 46 presents the resulting phasors for the closed-loop roll- 
spiral mode, ugg, which approximates the numerator "mode" u^. Using 
these phasors, vector polygons have been plotted for this mode (see 
Figure 46) from the equations of motion. Figure 13. Minor vectors 



Loop Closure 


Sstk 

< tk 

p/s 


s A 



nn 


LIT 



Figure 45. Closure of p/s ♦ $stk t0 Obtain Time Vectors, 
F-4J, a Q « 21 deg, - 5.5 deg 
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correspond to negligible terms; small polygons indicate little contri¬ 
bution to the modal response. Examination of these vector polygons 
leads to elimination of a number of terms from the equations of Fig¬ 
ure 13. Specifically, the kinematic equations are no longer needed, and 
so the matrix for approximation of the coupled numerator N t stk reduces 
to the form shown in Figure 47. The derivatives Y y , N$ aC k» and Z w are 
also minimal for this particular airframe and flight condition. How¬ 
ever, under other circumstances they can be quite large, so they have 
been retained to obtain a more general expression for 

The determinant of this matrix may be developed as a combination of 
uncoupled longitudinal and lateral-directional terms plus cross-coupling 
terms, as follows: 

<tk ; * «atk { + M q + M a) s ■ + Z wMq)] 



This is of the form 

where (^stk) unC0 upled con8 i 8ts of two second-order terras: a longitu¬ 
dinal term equivalent to the conventional short-period (denominator) 
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S - (Mq + Mi) 

-M 0 - MiZ w 



-Mg 


q 

-1 

8 - Z„ 





a 


-i-a 

^^stk 

- <£r 

-Xg 


P 


-NS, 


a - NJ 

-Ng 


r 



"“"I 

cos o 0 

s - Y v 


8 


Figure 47. Matrix for Approximation of Coupled Numerator 

approximation; and a lateral-directional term which is the 3 DOF N $ 8tk 
(with Y $ st k = 0) (Reference 20). The coupling "gain" terms are: 


and 


“ —MpjCg j"“"" cob a o 
6 stk 

K n - MgN£ cos a 0 


Servo analysis techniques can now be utilized to examine the influ¬ 
ence of each of the cross-coupling derivatives by writing the numerator 
in characteristic-equation form, i.e., 


stk uncoupled 


+ % 


and examining the migration of the roots as each of the "gains” is 
varied. 


Setting the coupling term K N to zero, the characteristic equation 
may be written as: 


1 






stk uncoupled 


0 
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Figure 48 is a system survey plot showing the effect of varying 
positive £ a on the uncoupled numerator The upper plot la a root 
locus and the lower plot a Bode-slggie (see Section II) showing the 
influence of increasing (negative) gain. These Indicate that the com¬ 
plex zero associated with the longitudinal short period, <^ p , i 8 driven 
to higher frequencies parallel to the jw-axie; and the conventional roll 
zeros, l/T^ and 1/T^ 2 , are driven apart. The influence of £ £ is found 
to be negligible in this case, but this Is due to a very small value of 
the ratio 6^ ( at the flight condition being evaluated) which 
is a factor in the "gain" term. 

The prime contributor to cross-coupling influence on the F-4J roll 
numerator is % ■ MgNg cos a 0 as shown in Figure 49, which is a survey 
for the equation 


1 



0 


where primes have been used to denote that K^> has already been ac¬ 
counted for, and double primeB indicate that both and K N have been 

accounted for. K N produces a locus opposite in direction to that of 
K/ , since both are pure gains but have opposite sign. The root loca¬ 
tion is much more sensitive to changes in Kjj, as the relative Bode 
amplitudes in Figures 48 and 49 show. 

The "gains" K/, and K N are functions of £&, N£, Mg, Ng^^, and 
*^stk* ^ile and have ^ een * traaaed * n the foregoing, it should 
be obvious that the influence of changes in magnitude of any of these 
derivatives can be equally represented by Figures 48 and 49. Changes 
in sign of any of the derivatives would reverse the direction of root 
migration in either case. 

The results shown here represent the usual derivative signs that 
prevail as AOA increases toward stall, l.e., £ £ positive, negative, 
and Mg negative. However, Mg depends upon several configuration details 
and can be either positive or negative (Reference 30). For the above 
signs, uncoordinated (8 j 4 0) flight results in reduced dutch roll 


7b 





(-.506H.943) [.386 ,.978| 



w or I<rl (rod/sec) 

inn 



e 48. System Survey of Effect of «£ ^ oi 
F-4j; - 0, a Q - 21 deg, 3 0 - 5.5 deg 












.gure 49. System Survey for Effect of on Ng , i 
£' a - 3.4251 1/sec 2 ), F-4J; a 0 - 21 deg, 3 0 “ S 5.5 








stability (Figure 43 or 44) and a corresponding shift of the roll 

numerator zeros. These shifts portend increasing closed-loop roll 

control problems with sideslip (Figure 45). 

The accuracy of the 5 DOF N $ 8tlt approximation is evident from 

Table 4: the coupled, 6 DOF numerator has been matched almost identi¬ 
cally. With this validation of the approximation, the following dimple 
coupled approximate factors were developed for this flight condition (it 
must be stressed that these literal approximate factors apply to this 

aircraft and flight condition only,and may not be accurate for other 
situations). For the general form 


< tt - («* + 2C,V + »$K. 2 + 2t{ p «{ p . + .* p ) 
: -<n; + Y v ) “ /NaMg cos a 0 

2 • * ... 

- Ng cos a 0 + cos a Q 

2 ^p^p " “CMq + Mi) + /N>p cos a Q 

4p : ~ M a + z w«q + 'N'M 3 cos a 0 


Table 4 includes the values given by these approximate factors, 
indicating that they also are accurate and therefore, encompass the key 
derivatives. 

These literal factors clearly indicate that longitudinal-lateral 
cross-coupling in is primarily through the product N^Mg. If 

either derivative is zero, the crosB-coupling terms disappear and the 
approximate factors become similar to the literal expressions for un¬ 
coupled lateral dutch roll and longitudinal short-period (Reference 20). 
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TABLE 4 


N ^stk NUMERAT0R APPROXIMATIONS 


NUMERATOR 


N 

; sp 

w ap 

6 DOF coupled 

-.531 

.909 

.778 

1.361 

5 DOF (approximate) 
coupled 

-.533 

.892 

.766 

1.370 

Based on literal 
approximate factors 

-.574 

.940 

.667 

1.544 


In addition, for non-zero Nj and Mg the aircraft coupling Is aensltlve 
to sideslip since 


3 2 N- . 

1 3g3a 


a* 

1 I z 


so that Figure 49 may be viewed as a locus of root migration with side¬ 
slip for constant C n g a . 

From the standpoint of vehicle design, any airframe changes which 
modify ilther of these cross-coupling derivatives can be expected to 
have a significant influence on high AOA flying qualities and departure 
characteristics. 

C. N? stab NUMERATOR 

Analysis of the influence of aerodynamic cross-coupling on the 
numerator for longitudinal control, N 6 sta ^» uses an approach identical 
to that taken with The first step is to make the development 

tractable by eliminating insignificant terms in the nine-by-nine matrix. 

It wao shown in Subsection II.A.3.b (Figure 19) that constraining 
velocity does little to change the zeros influenced by cross-coupling, 
i.e., I/Tq^, 1 /Tj|« and uij. Additionally, we find from Figure 50 that the 
equations for the Euler angles 4) and 9 can again be eliminated. This 
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Figure 50. Time Vectors for Closed-Loop Mode 1/T S 
(Approximating I/T 03 ), F-AJ; 
a - 21 deg, B - 5,5 deg 





still leaves a six-by-six matrix which is too large and complex to work 
with comfortably. 

Applications of time vectors, as was done with is accom¬ 

plished by closing a pitch rate loop with a feedback element of K/s 
(to obtain body-axis 0). Modal response coefficients are then calcu¬ 
lated for the closed-loop pole driven into 1/Tq 3> and the time vectors 
for each equation constructed as shown in Figure iO. Except for the 
a equation, few terms are small enough to be eliminated by inspection. 
Further,it seems unwise to restrict the general applicability of the 
results by eliminating the a equation. The main benefit of the time 
vectors is to demonstrate that the modal response for a root at 1/Tq 3 
contains considerably more lateral than longitudinal motion. 

Thus, it is necessary to expand the complete six-by-six matrix of 
Figure 51 in literal form, evaluate each of the terms, and eliminate 
products of small numbers. The resulting polynomial is then separated 
into uncoupled and coupled terras to form the approximation! 

«« Bt .b u ■ B w)|“ 4 ■ (*; + ■; + »„>■»+ |+ Y V > + HiYv + N0 dyn |*2 

- CO. cfc + /^Y v + Jj- X|3 dyn + (JfjNji - ■*» n o| B 

+ C0B “o + *>3 ■ ln °o|| f - <>»; + V s + co^\] 



where in the "dynamic axis" system of Reference 18, 

N Bdyn 2 N & cos “o " *0 8ln a o 


* 0 dyn 5 ^ B coe a o + N 0 sin a„ 
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Figure 51. Matrix for Approximation of Coupled 

(<fr 0 - Y 0 - 0, 0 - a, 8 0 Small) stab 


The terms not containing cross-coupling derivatives are seen to 
separate into a first-order and quartic function of a. The first- 
order Is the conventional literal approximation for the longitudinal 
zero 1/Tg 2 , i.e., -Z w . The quartic contains all lateral-directional 
derivatives and, with the exception of the sin a Q terms, is the same as 
the conventional body-axis uncoupled lateral characteristic equation 
(Reference 20). This polynomial thus factors into the parameters u^, 
1/T r , and 1/T S . 

The first cross-coupling term Involves as the "gain," a free s, 
and a quadratic containing only yawing derivatives. The second cross- 
coupling expression has N£ as a "gain,” a free s, and a quadratic con¬ 
taining only rolling derivatives. 

The kinematic terms Z p and Z r were negligible for N t 8tR » hut here 
dictate the influence of the aerodynamic crosB-coupling. As a result, 
coupling in the 0 * ^ 8 tab t,uraet,ator * 8 un explicit function of sideslip. 
It Is Interesting to note that Mg does not appear in the approximation; 
a quadratic coupling term involving M B z 6 8 tab^ VT o was found t0 he ne 8 1:t- 
gible. 
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The Influence of the Individual cross-coupling terms can again be 
observed by servo analysis methods, setting 






^stab ] 


and performing system surveys. The survey for 


f ^stab^ 


^ ^stab^ uncoup i ec j 


0 


is shown in Figure 52. The "poles" 1 /t| and 1 /Tq 2 couple to create a 
second-order closed-loop root, ug s , which is on the jw-axis. The dutch 
roll "pole" is moved to the left half-plane, while 1/T^ is moved slight¬ 
ly to the left. 

Figure 53 shows the additional movement when the cross-coupling due 
to Nq is included. The latter moves things further in the same direc¬ 
tion, However, it can be determined that the first-order coupled zero 

in the left half-plane comes from the roll subsidence mode and the zero 

farthest into the right half-plane (previously called l/T^) comes from 
a coupling of the spiral mode and the. conventionale 1 /Tq 2 • Since the 
time vectors indicated this mode is dominated by lateral-directional 
motion it is now labeled U/Tg)". The prime notation in Figures 52 
and 53 is used to keep track of the number of "loops" being closed 

and the influence of each "closure." From here on the prime notation 
will be dropped and the lateral-longitudinal coupled zeros identified 
as t*>®, 1 /Tjj, 1 /Tg, etc. 

Since each of the cross-coupling terms has similar Influence on 

the migration of zeros, no insight is lost by combining the two coupling 
terras via the concept of the previously noted "dynamic" terra Then 
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The system survey for the combined coupling influence on ^ 6 ata ^ is shown 
in Figure 54; the coupled roots, of course, are the same. 

Finally, it should be noted that aerodynamic cross-coupling in the 
numerator is a direct function of sideslip squared . Since 3C„ is 
proportional to Cj^ a (e 0 ) and is proportional to C^CSo), the "gain" 
term is 


*a dyn 0o ~ (% a e 0 c °s a 0 + OS^o sin a 0 )B 
~ cos °° + Sa sin a °) e o 

Thus the zero 1 /t| is drastically influenced by static aerodynamic 
cross-coupling and magnitude of sideslip. 

D. SUMMARY 


Literal analytic expressions have been developed to show the influ¬ 
ence of static aerodynamic cross-coupling derivatives and asymmetric 
flight on denominator and numerator roots. The analysis has shown that 
the shift in denominator roots depends upon all three static cross¬ 
coupling derivatives, NJ, and M 3 . The direction of root migration 

is controlled by the sign and magnitude of MpN^/ «C^. The magnitude of 
root migration is determined by ^0" 
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For the F-4J the. principal cross-coupling Influence on the roll 
numerator, N t 8tk » ia due to the term cos a 0 . The literal approxi¬ 

mate factors which describe the coupled numerator roots are 

“(Nr + Y v > - cos n 0 

Ng cos a Q + cos a Q 

“(Mq + M$) + /N£Mg cos a 0 

-M a + ZyMq + /K'M $ cos~a 0 

The pitch numerator, N 6 8ta jj» ls influenced by a more complicated 
combination of kinematic and aerodynamic terms. Existence of a first- 
order RHP zero depends upon 


- 


4sp u sp 


/a N 0 “ ^B n o “ Y v(^r cos °o " sin a o) 
Jt' cos a 0 + N' sin a 0 


The above approximations may not apply to configurations significantly 
different from the F-4J. 
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SECTION IV 


CONFIGURATIONS SELECTED FOR SIMULATION 


The seven aerodynamic derivatives which have been shown to dominate 
open- and closed-loop dynamic characteristics at high AOA are summarized 
in Table S. These are comprised of six static stability and cross¬ 
coupling derivatives and one dynamic damping derivative. Together they 
limit the useful flight envelope via roll reversal and wing rock tenden¬ 
cies, and they limit the safe flight envelope via departure suscepti¬ 
bility. One goal of the piloted simulation was to demonstrate this 
premise. 


TABLE 5 

KEY MANEUVER-LIMITING PARAMETERS 




Negative N$ a or N$ D 

Roll reversal 

Positive M 0 

Pitch up 

£ %» £ at £■ p | 

( Wing rock 

Ng. K I 

\ Nose slice 


V Roll divergence 


This section presents the vehicle aerodynamic configurations se¬ 
lected for the piloted simulation and the predicted high AOA dynamic 
characteristics. 

Four "configurations” were obtained by altering the key aerodynamic 
coefficients of a baseline. The three altered configurations incor¬ 
porated changes in, respectively, Cjj p (a), C^(a,fi), and C n (a,&) combined 
with a change in C m g, to provide specific alteration of key open- and 
closed-loop parameters discussed in the preceding two sections. The 
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attempt was to influence the denominator only, the numerator only, and 
the combined denominator and numerator, thereby to provide distinct 
variations in high AOA departure susceptibility, onset warning, and 
motion severity. The influence of artificially altering aerodynamic 
coefficients via a lateral-directional augmentation system was also 
investigated with two of the configurations. 

A. AERODYNAMIC CONFIGURATIONS 

1. Configuration A 


The baseline case is ider ‘led as Configuration A. This F-4J- 
based configuration has been thoroughly reviewed in the preceding sec¬ 
tions. The key open- and cloaed-loop dynamic parameters are summarized 
In Figure 55 for comparison with the other configurations discussed 
subsequently. Dutch roll and numerator root migrations are shown for 
increasing a at zero 5 (solid lines) and at 5.5 deg B (dashed lines). 
The solid triangle (A) indicates the AOA above which C n ^ becomes 
negative at zero 0. Recalling the discussion in Section II about the 
shortcomings of linearized frozen point analysis, the important aspects 
to note here are: 

e The low dutch roll damping at AOA below C n _ - 0 
(wing rock). ^yn 

• A or 1 /Tj, negative at AOA above 18 deg or at 
small 0 (roll reversal). 

e 1/T*, and 1 / T s both large and negative when S i* 0 
(nose slice parameters). 


Dutch roll and roll numerator roots sensitive to 
sideslip. 


The plots in the lower left corner indicate the magnitude of the 
aerodynamic coupling terms evaluated at 5.5 deg B. In the AOA range 
from 20 to 30 deg, jCa is of approximately the same magnitude as but 
is of opposite sign. In this same AOA range, Is of about the same 
magnitude as Njj and of the same sign. For this relatively mild sideslip 
condition (5.5 deg), the static cross-coupling derivatives are of equal 
magnitude to the basic static derivatives. 
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2. Configuration B 


For Configuration B the roll rate damping coefficient, C£ p , was made 
large, negative and essentially Invariant over the range 10 < a < 35 deg 
(see c ip 2 » Figure 56). The object of this change Is to Increase the 
dutch roll damping and hence reduce the wing rock tendency. Since all 
other parameters remain unchanged, this should not alter the basic roll 
reversal tendency or the nose slice departure characteristics; but It 
should eliminate or reduce wing rock as a warning of impending depar¬ 
ture, However, note again in Figure 57 the sensitivity of the dutch 
roll mode to sideslip. At B » 5.5 des there is little difference 
between Configuratlona A (Figure 55) and B (Figure 57), which indicates 
that the influence of «Cp is small with respect to that of jCjj. 



Figure 56. Change in for Configuration B 
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3. Configuration C 


For Configuration C the roll due to sideslip coefficient, Cg^, was 
altered (see Figure 58) to approximate that of the F-14A over the range 
15 < a < 45 deg, taking into account the difference in roll inertias of 
the two aircraft, i.e., 




C *Sp-U 


r *F-4 

X *P-U 


Figure 59 presents the parameter plots for Configuration C. The 
upper left figure shows that the change introduced in Cj^g and Cg a haB 
shifted the dutch roll mode toward the left. It remains in the left 
half-plane for 8 ■ 0 and only bulges into the right half-plane over a 
small a range at 0 ■ 5.5 deg. A major difference in thiB configuration 
occurs in the roll subsidence and spiral modes, which in the previous 



Figure 58. Change in C» for Configuration C 

p 
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cases have remained on the real axis and therefore have not been shorn. 
In this case, 1/T R and 1/T 8 couple to form a "lateral phugoid” mode, 
u, tRI as AOA increases. At 36 deg AOA the frequency of this mode is 
approximately twice that of the dutch roll, with approximately the same 
damping. Comparison of Figures 35, 39, and 59 show this to be a valid 
representation of F-14A root locations, but mode identity may be differ¬ 
ent. 

Comparison indicates that increasing 6 shifts the u>g R locus toward 
the right half-plane. However, both the dutch roll and lateral phugoid 
modes are now relatively unaffected by 3 . 

The upper right-hand plot for is essentially the same for Con¬ 
figuration C as that for Configurations A and B. While this configura¬ 
tion has better open-loop stability, the closed-loop characteristics 
should be little different from the previous two configurations because 
the roll numerator roots are essentially unchanged. 

Absence of the solid triangle symbol indicates that ia P 0Bi “ 

tive throughout the AOA range. 

The lower left-hand figure shows that the term (which derives 
from the variation of Cj^ with a) is small and negative. Since this 
configuration also has a large increase in £jj (see Figure 58), the 
aerodynamic cross-coupling derivative *Cq 1b now muc h waller than | 
throughout the angle-of-attack range. Comparing the values of £ £ for 
Configurations A and C with the differences in denominator root migra¬ 
tions with B makes it apparent that is the dominant purameter in 
shifting the Configuration A dutch roll mode into the right half-plane 
when B 0. The lower right-hand plot of Figure 59 indicates that the 
change in ttnd ^8 ha ® decreased 1/T§ at AOA above 20 deg. 
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4. Configuration D 


For the fourth case the yaw due to sideslip coefficient, C n ^, was 
modified to approximate that of the F-14 as shown by the dashed line In 
Figure 60. In addition, for this case only, the sign of C m g was re- 
versed to provide positive pitch with sideslip. As predicted by the 
approximate factors of Section III, the result Is an airframe quite 
Insensitive to sideslip (Figure 61) and, most pertinently, having numer¬ 
ator roots (l/T^ and 1/T§) which do not penetrate as far into the right 
half-plane. 


Roll reversal is delayed until 
tive until a * 25 deg. 


> 20 deg and C, 


n ®dyn 


remains posi- 



Figure 60. Change in C n g for Configuration D 
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5. Anticipated High AOA Maaeuver- 

Llnltlng Characteristics 

Review of the preceding four configurations should go a long way 
toward sorting out the relative Importance of denominator (open-loop) 
vs. numerator (closed-loop) characteristics on high AOA departure sus¬ 
ceptibility and severity. Configurations A, B, and C have nearly iden¬ 
tical RHP numerator root locations and therefore should have similar 
closed-loop divergence tendencies, although the motion (l.e., roll vs. 
yaw) may differ somewhat due to Cj^, C^g, and Cj, a changes. On the other 
hand, the Configuration D numerator roots do not penetrate so far Into 
the RHP and are relatively insensitive to sideslip. Thus one might 
expect this vehicle to appear less susceptible to or have less severe 
departure. 

Alternately, if departure susceptibility or tendency should turn out 
to be strongly influenced by dutch roll (open-loop) root location, then 
one would expect Configurations C and D to be similar since both are 
near neutral stability and relatively insensitive to sideslip. Any 
differences between Configurations A and B would be dependent upon side¬ 
slip. 

The dynamic stability parameter, C n _ , has been referred to pre- 
R dyn 

viously. The AOA at which this parameter becomes negative was identi¬ 
fied in Figures 55, 57, and 61. Figure 62 presents plots of 
C n g^ vs. a for the four configurations. Based upon the criterion 

of C na > 0.002 for departure resistance (Reference 31), then only 
Pdyn 

Configuration C should be resistant, Configurations A and B should be 
quite departure prone, and Configuration D might be mildly departure 
prone. 

Another criterion, Reference 32, employs both open- and closed-loop 
parameters, l.e., 
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Figure 62. 


C no for Configurations Simulated 
s dyn 


For a number of years It has been recognized (Reference 33) that unde¬ 
sirable ratios lead to roll reversal and pilot-aggravated wing 
rock (PIO). However, in Reference 32 Weissman plots c ng dyn vfl ’ LCDP and 
identifies four regions of increasing departure and spin susceptibility 
and severity (Figure 63) based upon experience with a number of airframe 
configurations. The loci of our unaugmented airframe parameter values 
over the range of 16 < a < 35 deg are shown in Figure 63; circles repre¬ 
sent Configurations A and B (since C* p has no influence on either param¬ 
eter), squares represent Configuration C, and triangles Configuration 
D. The criteria predict high departure/spin susceptibility with strong 
rolling departures for Configurations A and B, low spin susceptibility 
and mild rolling departures for Configuration C, and no departures for 
Configuration D. Thus, all regions of departure susceptiblity and 
severity are exercised with the aerodynamics selected. 

The Figure 63 prediction that Configuration D should be less suscep¬ 
tible to departure is in agreement with the first assessment above based 
on numerator root locations. 
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B. AUGMENTED CONFIGURATIONS 


Lateral-directional augmentation was synthesized for Configuration A 

(the most departure prone) and for Configuration C (C n always posi- 
. ^dyn 

tive). The augmented configurations are identified as and C 2 , re¬ 
spectively, while the unaugmented configurations will be referred to as 
and Cj. The augmentation mechanization is representative of that 
used in current high-performance fighters to improve both open- and 
closed-loop handling parameters, i.e., washed-out yaw rate and lateral 
acceleration feedbacks to rudder to increase dutch roll damping; roll 
rate feedback to differential stabilator to reduce or eliminate wing 
rock and augment dutch roll damping; and lateral stick to rudder cross- 
feed (SRI) to reduce or eliminate adverse aileron yaw for pilot inputs, 
thereby providing a positive shift in the closed-loop parameters LCDP, 

«J, ^ 1/V 


Undesirably, the roll damper also produces adverse yaw when opposing 
uncomraanded roll; it may increase nose slice departure tendencies and 
pro-spin surface deflections. A block diagram for the augmentation 
mechanization is shown in Section V. 

The root migration surveys for Configuration A 2 are shown in Fig¬ 
ure 64. Comparison with Figure 55 shows that at 0 - 0 the dutch roll 
damping Is increased considerably; it does not go unstable until AOA 
exceeds that for C n ^ becoming negative (A symbol). Also, the AOA at 
which ( 1 )^ becomes negative has heen increased approximately 5 dt g and the 
negative values of both 1/T^ and 1/T® have been reduced. Thus, within 
the augmentation system authority limits and near zero 0, Configura¬ 
tion A 2 stability is increased and departure susceptibility decreased up 
to about 23 deg AOA. However, for 0 0 the augmentation benefits are 

overcome and the dynamic parameters revert back to air ,: r«irae-alone 
values. 
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Figure 64. Key Parameters for Configuration A 2 


Root migration survey plots for Configuration C 2 are presented in 
Figure 65. Comparison with Figure 57 again shows Improvement in dutch 
roll stability. However, the most dramatic change is in u>^ which now 
remains positive up to 35 deg AOA (but note is small and negative 

beyond 31 deg AOA). Thus roll reversal problems are essentially elimin¬ 
ated and numerator root sensitivity to sideslip is decreased consider¬ 
ably — again subject to augmentation system authority limits. 

The change in departure and spin characteristics cannot be pre¬ 
dicted for these configurations using the static aerodynamic parameters 
of Figure 63 because the augmentation system influence is frequency 
dependent. But a somewhat similar plot can be constructed using the 
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Figure 65. Key Parameters for Configuration C 2 
2 * 

parameters and Such a P lot *•» presented in Figure 66 for 

both unaugmented and augmented versions of Configurations A and C. The 
2 

Intercept of the w<j> axis with Cd^d has been arbitrarily shifted to en¬ 
hance similarity for the unaugmented configurations with Figure 63. 
Note that the root plots for the unaugmented Configurations Aj and 
(solid lines) are quite similar to those of Figure 63. The plots 

(Figure 66) for Configurations A 2 and C 2 (dashed lines) are moved fur¬ 
ther into the upper right quadrant, which is identified as departure 
resistant in Figure 63. Similar interpretation would forecast A 2 to 
be departure resistant up to about 24 deg AOA but to have moderate to 
severe departure tendencies at higher AOA, while C 2 should have no 
departure tendencies. 


w n and u)^, it is move 
0 dvn 

appropriate to plot because the changes in sign (or lack thereof) 


Despite the direct relationship between C n; 

C d w d 

i are in closer agreement. 


and { 


p dyn 
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Figure 66. Versus C d u> d for Augmented and Unaugmented 
Configurations A and C 


C. SUMMARY OF CONFIGURATIONS AND 

PREDICTED HIGH a CHARACTERISTICS 

Table 6 summarizes the six-configuratlon matrix employed in the 
piloted simulation. For A^ a sequence of roll reversal, wing rock, nose 
slice, and finally strong rolling departure with high spin suscepti¬ 
bility is predicted with increasing AOA. For A 2 the augmentation system 
and crossfeed are expected to minimize or eliminate the roll reversal 
and wing rock predeparture warnings. The predicted characteristics are 
nose slice followed by moderate rolling departure. For B, predicted 
characteristics are roll reversal, nose slice, and strong rolling depar¬ 
ture. Warning In the form of wing rock should not be present due to the 
large roll damping of this configuration. Note that Configurations A 2 
and B allow comparison of high AOA stall/departure characteristics with 
high aerodynamic roll damping (Configuration B) and with artificially 
augmented roll damping (Configuration A 2 ). 

The unagumented Configuration Cj is predicted to exhibit roll rever¬ 
sal, wing rock, and mild rolling divergence with increasing AOA but to 
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TABLE 6. CONFIGURATION MATRIX 


CONFIGU¬ 

RATION 

FCS 

AERODYNAMIC VARIANT 

PREDICTED 

CHARACTERISTIC 

A 1 

Basic 

Basic F-4J 

Wing rock 

Roll reversal 

Nose slice 

Roil departure 

A 2 

Augmented 

Nose slice 

Roll departure 

B 

Basic 

Increased Co 

P 

Roll reversal 

Nose slice 

Roll departure 

C 1 

Basic 

Increased Co_ 1 

B [ 15 < a < 45 
Decreased ) 

Wing rock 

Roll reversal 

Roll departure 

c 2 

Augmented 

Cna, > 0 
None Pd y n 

LCDP > 0 

D 

Basic 

Increased C n . ) 

6 > a > 15 
Decreased C na ) 

Positive C ffl g 

Wing rock 

Roll reversal 

Pitch up 


have no nose slice. Augmentation is employed with C 2 to improve maneu¬ 
vering control and to determine if it will degrade the departure resis¬ 
tance of this configuration. On the basis of C na and (or LCDP), 

Pdyn v 

no departure tendency should be anticipated for C 2 . Finally, D should 
exhibit roll reversal and wing rock warnings but little lateral-direc¬ 
tional departure tendency. However, the change in sign of to pro¬ 
vide positive pitching moment with sideslip should result in pitch-up, 
which would be expected to aggravate any high AOA departure tendency. 
All other configurations have negative C^, should pitch down with 
increasing sideslip, and hence should require more effort of the pilot 
to maintain high AOA (longitudinal stick cue). 










PILOTED SIMULATION 


The purposes of the piloted simulation were to validate the influ¬ 
ences of the key aerodynamic coefficients in determining departure 
characteristics, evaluate the influences of varied maneuver limiting 
factors on high AOA maneuvering control, and identify potential flying 
qualities criteria. While aircraft motion is highly desirable for such 
validation and assessment, departure and post-stall gyration (PSG) 
severity obtained with the 6 DOF analytic models indicated that a 
moving-base simulation would offer very little benefit. Noee slice 
motion was so rapid that 130 dog/sec yaw rate could be reached within 
two seconds after departure onset; peak angular accelerations as large 
as 

p ■ 500 deg/sec 2 

q - 200 deg/sec 2 

r - 150 deg/sac 2 

were obtained. Washout requirements to prevent hitting displacement 
y stops would have to be so rapid as to negate motion benefits. There- 
! ( fore, a fixed-base simulation was selected having an unrestricted out- 
of-cockpit horizon display capable of such rates and accelerations. The 
aerodynamic model was also expanded to accommodate inverted flight 
maneuvers. 

This section presents an overview of the simulator, Including dis¬ 
plays, airframe, and flight control system modeling, cockpit, and force 
and moment equations. 
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A. SIMULATOR 


The simulation was performed in the McDonnell Aircraft Company's 
fixed-base Manned Air Combat Simulator dome, a 20 ft hemisphere Identi¬ 
fied as MACS-1. Physical aspects of the simulation are summarized in 
Figure 67. The horizon and target are projected on the inside of the 
hemisphere. The-, cockpit is located at the center of the dome. The out- 
the-window, head-up (HUD) and head-down (HDD) displays, and cockpit 
layou: are sb indicated in Figure 67. Seat cues consisted of normal 
acceleration and buffet motion provided through an inflatable seat 
bladder. The frequency and amplitude of the buffet oscillation varied 
as the aircraft angle of attack increased, starting at about 14 deg. 

A TV, projection of a gimballed model target aircraft provided a 
maneuvering tracking task. Target motion equations provided realistic 
maneuvers while allowing the aircraft to fly well beyond the departure 
region for the test aircraft. The target was controlled by the computer 
operator via a special control panel. 

B. HEAD-UP DISPLAY 

The head-up display (HUD) provided information on aircraft attitude, 
heading, airspeed, altitude, and trim conditions. Figure 68 illustrates 
the HUD presentation as originally installed. Display of normal accel¬ 
eration provided a reference for trimming longitudinal control to 1 g 
prior to the start of a run. Trim thrust was set by adjusting throttles 
until change from trim thrust registered near zero. The Flight Path 
Angle Ladder displayed both flight path elevation angle and roll atti¬ 
tude. Orientation of the total velocity vector was provided through 
the Velocity Vector Indicator; when the velocity vector was outside the 
HUD field of view, the Indicator remained at the edge of the HUD and 
blinked. 

At pilot request, the HUD was simplified late in the simulation 
program. Airspeed, Heading, and Altitude bar displays were removed. 

The fixed reticle tracking gunsight was a part of the HUD. 
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MCDONNELL. AI ft C ft A 


6 Real Horton and Missile Projector 

7 Real Target Focus Lenses 

8 Real Targut Mirrors 

9 Real Target Projector 

10 Sound and Electronic Equipment 

11 Pit Area 


m-i m 


Displays: 

Horizon - 360 deg 4>, 0, ij) 

HUD - CAS, h, velocity vector 

HDD - $, 8, <|>, <*, M, etc. 

Sight - fixed reticle 

Cockpit: 

Basic F-4 

Seat Cues: 

Load factor, buffet 

Target: 

Gimballed model TV projection 

Figure 67. 

Manned Air Combat Simulator I 






Aircraft Heading (deg) 



Thrust Lever 

Normal - *-NZ 0.9 

Acceleration 

<fl) 


Figure 68. Heed-Up Display (HUD) Used in Simulation 


C. FLIGHT OOHTML SYSTEM 


The manual flight control mechanization (Btick, pedals, surface de¬ 
flections and rates, etc.) represented the basic F-A aircraft (e.g., see 
Figure 2). 

A special lateral-directional stability and command augmentation 
system (SCAS) was mechanized as shown in Figure 69. The SCAS feedback 
gains are production F-A values. The stick-to-rudder crossfeed gain 
and shaping were optimized for 10 < a < 20 deg and faded to zero at 
a < 10 deg to prevent adversely influencing the low AOA handling quali¬ 
ties. 


D. EQUATIONS OF MOTION 

Slightly modified versions of the standard McDonnell Aircraft Com¬ 
pany piloted simulation equations of motion were used. These are essen¬ 
tially the same as those presented in Appendix I, Part III, i.e., moment 
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equations in body centerline axes, force equations In wind axes, and 
aircraft orientation in standard Euler axes. The only change required 
was to expand the small angle a and & assumptions and to eliminate 
the tan 3 and cos - *’ 0 expressions of the d equation which gave discon¬ 
tinuity problems as 3 approached 90 deg. This required substituting 
the drequation for 4, i.e., 


4 - q - tan 3(p cos a + r sin a) + Zy/CmVip cos 3) 

■ qV T cos 3 - V T sin 3(p cos a + r sin a) + 2^/m 


and using the simple expressions a ■ sin”* 1 (w/Vj) and 3 ■ sin”*’ (v/V T ) 
for entry.into the aerodynamic look-up tables. 

E. AERODYNAMIC MODELS 

Since the major interest in the simulation was the analysis of 
departure onset, particular attention was given to realistically model¬ 
ing aerodynamic data for the range of 0 < a < 45 deg and 3 < ±30 deg. 
The basic aerodynamic coefficients were selected with the goal of pro¬ 
viding a reasonably accurate dynamic model of the F-4J over this o and 3 
range. To allow greater aircraft motion freedom typical of severe 
departures, post-stall gyrations, and spins without encountering unreal¬ 
istic discontinuities in the data, certain approximations were made for 
a > 45 deg, a < 0 deg, and 3 > ±30 deg, i.e., 

• For a > 45 deg the coefficients were faired to 
110 deg, then extended to 180 deg by either main¬ 
taining the value at 110 deg or fading the coef¬ 
ficient to zero at a - 180 deg. 

• The coefficients were assumed to be similar in 
value for positive and negative A0A so that all 
derivatives were taken to be even functions of A0A 
(i.e., mirror images about a - 0) and the coeffi¬ 
cients and C m to be odd functions of AOA. 

• Variation of the coefficients with 3 was Assumed 
to be linear for 3 > 30 deg. 

The resulting aerodynamic, coefficients are shown in Figure 70. 
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It is recognized that the assumptions made above could result in 
considerable error for extreme angles of attack and sideslip; however, 
the use of full-ranged data allowed the simulation to be operated 
throughout a violent departure without demanding an immediate return to 
initial conditions when lower aerodynamic data limits were reached. 

The aerodynamic data were stored in the digital computer as look-up 
tables. Total forces and moments were computed digitally and fed to the 
analog systems of the simulator to drive the cockpit displays. 

Physical dimensions of the airplane were based upon the P-4J. Iner¬ 
tias and weights are summarized in Section IX. 

F. DATA RECORDED 

Parametric data were recorded via three eight-channel Brush record¬ 
ers. All Euler angles, body axis rates, and control deflections were 
recorded, as well as accelerations, altitude, velocity, and thrust set¬ 
tings. 

For a limited number of runs the total body axis accelerations 
(p» q* 0 were recorded and inertial and aerodynamic components of 4 
were also recorded separately. 

G. TASKS AND MANEUVERS 

Piloting tasks were selected which exercised the open- and closed- 
loop departure parameters of interest. These were divid«i into two 
phases. The first consisted of familiarization maneuvers typical of 
"feeling out" stall/departure motions, warning, etc., of a new airframe. 
These are modeled after the suggested test matrix of Reference 34 plus 
an additional aggravated input maneuver for determining departure sus¬ 
ceptibility, proposed in Reference 35. The second phase consisted ot 
tracking tasks and maneuvers as suggested in Reference 36. 





1. Handling Fuillcriutlon and Aaaesment 

Five basic familiarization (F) maneuvers were employed: 

• Straight-ahead stall (Fj): 1 g , fixed thrust, 

longitudinal stick pull to produce d * 1 deg/sec 
until stall departure. 

• Constant attitude stall (Fo) J 1 g. wings level, 
holding constant pitch attitude while slowly re¬ 
ducing thrust until stall/departure. 

• dank-to-bauk turns (F 3 ): Constant altitude, 

±60 deg bank-to-bank turns at increasing AOA until 
departure; rudder coordination optional. 

• Wind-up turn (F^)s Constant thrust, wind-up turn 
to stall departure. 

• Full-stick-deflection maneuver (Fe)r From 60 deg 
bank, rapidly ramp full aft longitudinal Btick 
followed by full lateral stick out of the bank. 

After 8 sec neutralize controls. 

With the exception of maneuver F 5 there was no restriction on the 
type or magnitude of controls used by the pilots. In fact, pilots were 
advised to investigate various control techniques juBt as they would in 
a stall evaluation of a new aircraft. Maneuver F 5 (from Reference 35) 
was evaluated as an indicator of departure susceptibility. 

2. Precision Tracking Evaluation 

The target aircraft <a TV picture of a model airplane, projected 
onto the dome at the proper aspect and size) was programmed to perform 
three basic maneuvers during a slow pull-up. Target climb rate was 
fixed to produce 4*1 deg/sec to guide the subject pilots into repeat¬ 
able tracking runs of approximately 30 sec minimum duration before 
departure. The tasks started with the subject aircraft trimmed at 
u - 12 dig, essentially in trail. The pilot attempted to stay in the 
target's 6 o'clock position as if trying to overextend a firing 
solution. The tracking (T) maneuvers were: 
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• Straight-ahead pullup (Tj): Conatant-thrust, 

straight-ahead pullup keeping pipper on target 
until aircraft no longer controllable. 

• Climbing roll-reversals (T 2 ): Track target air¬ 
craft through aeries of climbing roll rever¬ 
sals ($ * ±40 deg) until aircraft no longer con¬ 
trollable . Constant'thrust. 

• Wind-up turn (Tj): Starting wings-level, track 
target into con«tant-thrust 60 deg bank wind-up 
turn until airciair no longer controllable. 


3. Departure Recovery 

Recovery techniques employed were at the discretion of the pilots, 
based upon their experiences in fighter aircraft. Recoveries were ini¬ 
tiated either: 

• Immediately upon positive Indication of uncora- 
manded motion (Rj), or 

• After about a 3 aec delay (R 2 ). 

Recoveries were also attempted in the SCAS-on configurations (Cases A 2 
and C 2 ) by turning the SCAS off upon departure (R 3 ). 

H. FLYING QUALITY RATING SCALES 


Previous etall/departure simulations for flying quality assessment 
(e.g., Reference 5) have shown the familiar Cooper-Harper scale (Fig- 
.--e 71) to be Inadequate because departure is centered on the first step 
i v the decision tree, i.e., controllability Instead of performance. It 
was preordained for this simulation that control would be lost (CH ■* 10) 
because the configurations were specifically selected to provide differ¬ 
ing departure warning, severity, and recovery In order to Identify gra¬ 
dations in these attributes. Therefore, one requirement of this program 
was to develop and validate a more appropriate stall/departure flying 
quality rating scale. 

The Reference 37 experimental program succeeded in developing rating 
scales for a similar loss-of-control and recovery situation Induced 
by flying into a vary strong wake vortex. This experience indicated 




Highly desirable desired perfoi 
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that each attribute must be evaluated individually and then the overall 
hazard assessed. Figure 72 is a tentative departure rating scale adapt¬ 
ed from the wake vortex experience. Here the stall/departure attributes 
are divided into four open-ended scales. The first is departure onset 
warning — clarity or lack of warning. Next, the attitude and attitude 
rate excursions are evaluated. Recovery characteristics are separately 
evaluated (whether the recovery was rapid and achieved by simple control 
application or whether it degraded down to the point of no recovery). 
Finallyi an overall assessment of hazard is rated. 

In anticipation that overall hazard assessment (acceptability) might 
vary with aircraft primary mission, provision is made for indicating 
scalar differences between training and operational squadron usage. For 
example, an operational aircraft might have very clear warning of depar¬ 
ture onset and extreme excursions if the pilot persisted, and be diffi¬ 
cult to recover, but be given a relatively good hazard rating, i.e., a 
2, on the basis of the very clear and consistent warning. Similarly, a 
training aircraft might have little warning of departure onset and rela¬ 
tively large attitude excursions, but recover very easily and quickly by 
merely releasing the controls, and again might have a relatively good 
overall hazard rating, such as a 2, on the basis of recovery. 

As an adjunct to the rating scales, the pilot is also asked to pro¬ 
vide qualitative information as to the nature of the aircraft motions 
observed. 

Although the rating scale of Figure 72 was conceived to cover both 
open- and closed-loop situations and specifically is not spin oriented, 
it was employed in a separate departure/spin susceptibility piloted 
simulation (Reference 38) to obtain an early assessment of its useful¬ 
ness. Results were somewhat inconclusive since the Reference 38 task 
involved preprogrammed, open-loop pilot inputs to assess the departure/ 
apln susceptibility of an airframe having variations in surface author¬ 
ity and actuator rate limiting characteristics. Nevertheless, valuable 
comments were obtained (Reference 39) along with a recommended expansion 
and restructuring of the pilots' qualitative assessment to provide 
greater flexibility and provoke additional commentary to identify what 
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DEPARTURE ATTITUDE AND ATTITUDE 



• Wfiat was primary departure 
variable? 








the pilot is rating. Accordingly, a losa-of-control/departure/recovery 
debriefing guide was prepared, Table 7, to augment the rating chart. 
This guide addresses the four topics of the Figure 72 rating scale plus 
an additional section devoted to post stall gyrations (PSG) characteris¬ 
tics and severity* Provision is made under each topic for the numerical 
rating (STI) from the Figure 72 scale. The pilot is requested to answer 
those questions appropriate for the maneuver flown. Most questions per¬ 
tain to characteristics of the vehicle motion and are constructed to 
encourage simple yes/no, multiple choice, or short written answers. 
Cooper-Harper (CH) ratings are also requested for the two tasks which 
Involve pilot effort to exert control on the vehicle (l.e., prevent 
departure and accumpllsh recovery). Finally, it should be noted that 
the questionnaire was prepared for use in flight as well as simulation. 
Thus, some questions are not appropriate for fixed-base simulation 
(e.g., assessment of vehicle acceleration). 






TABLE 7. LOSS OF CONTROL/DEPARTURE/RECOVERY 
DEBRIEFING GUIDE 1 . 





















SECTION VI 


HIGH AOA FLYING QUALITIES ASSESSMENTS 


Effects of the maneuver-limiting factors on high-angle-of-attack 
flying qualities were assessed by three pilots utilizing the prescribed 
set of non-tracking and tracking maneuvers. The two principal pilots 
ware highly qualified U.S. Air Force fighter test pilots, graduates of 
the Air Force Test Pilot School at Edwards AFB, with a broad range of 
practical experience. The first pilot, RC, has flown a varied selection 
of aircraft, from the F-4 to the B—52; the second pilot, JF, was experi¬ 
enced in operational fighter and attack aircraft, including the A-9 
stall/spin evaluation, and la a former instructor at the Test Pilot 
School responsible for stAll/spin curricula. The third pilot, RH, was a 
handling quality engineer/general aviation pilot who also accomplished 
all of the Initial simulation debugging, procedural checkout, primary 
pilot briefingB and debriefings, etc. 

The combination of airframe configurations, familiarization and 
tracking maneuvers, and recoveries resulted in the 98-cell matrix of 
Table 8. This matrix was flown first as a sequence of runs progressing 
horizontally across the various maneuvers for a given airframe. After 
the complete matrix had been accomplished once, a second scries was per¬ 
formed progressing vertically through the airframe configurations for 
each tracking task. 

A total of 1088 formal runo were conducted. These were fairly 
evenly divided between the various configurations and pilcfcs, aa indi¬ 
cated by Table 9. 

One of the first tasks of the simulation was to validate or modify 
the ratings scale and debriefing questionnaire to the satisfaction of 
the evaluation pilots. The final versions were used throughout the 
remainder of the simulation. Accordingly, this section starts with the 
rating ocale assessment and then presents the airframe configuration 
assessments. 
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TABLE 9. NUMBER OF RUNS 


CONFIGU¬ 

RATION 

PILOT j 

TOTAL 

RC 

JF 

RH 

A 1 

59 

79 

79 

217 

a 2 

68 

48 

79 

195 

B 

61 

94 

58 

213 

C l 

44 

77 

40 

161 

c 2 

53 

62 

58 

173 

D 

35 

62 

32 

129 

TOTAL 

320 

422 

346 

1088 


A. RATINGS SCALE AND QUESTIONNAIRE ASSESSMENT 

The debriefing questionnaire was answered at the completion of the 
familiarization runs and again at the completion of the first set of 
tracking runs. During this period modifications were worked out for the 
rating scale. The modified scale was then used for the final set of 
tracking runs. 

1. Questionnaire 

The debriefing questionnaire proved essential to obtaining a consis¬ 
tent set of information for interpretation of results. However, the 
simple yes or no answers requested did not suffice. The tasks and 
resulting aircraft motions were so complex as to require considerable 
explanation. For example, more than a single mode of departure or spin 
was observed in some cases. Departure warning and severity also varied 
with rate of stall onset (&) and pilot usage of controls. Thus accom¬ 
plishment of the questionnaire often became quite time-consuming, with 
several pages of accompanying explanation. This was augmented by tape 
recordings made, on-line, of pilot commentary while performing the simu¬ 
lation runs. As a consequence, pilot responses were so lengthy — and 
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in so many cases repetitive - that it is not practical to incorporate 

them in this report. 

However, the basic conclusion reached by both pilots and experi¬ 
menters was that such a questionnaire is highly desirable for stall/ 
departure/recovery investigation. The only change in format deemed 
necessary is provision of more space to accommodate lengthy answers. 

2. Rating Scales 

As in the previous simulation (Reference 5), the pilots felt that 
the Cooper-Harper rating scale was inappropriate for stall/depatture/ 
recovery evaluation where aircraft characteristics rather than task 
performance are being assessed. On the other hand, the overall concept 
of the Figure 72 departure ratings scale was appreciated, although som^ 
alterations were required before this scale met with full approval. In 
particular, the pilots did not like the open-ended aspects of the scales 
(the scale extends beyond the first and last descriptors). This clearly 
gave problems, and they preferred that the scale be bounded with a 
simple descriptor at each end. Second, the time required to recover and 
the control application necessary to effect recovery needed to be sepa¬ 
rated. Recovery controls also deserved a further breakdown, with one 
rating for required control complexity and a second rating for control 
application timing. Finally, it was decided that an assessment of 
hazard was difficult to make because this changes with situations such 
as low versus high altitude. The pilots preferred to provide an overall 
assessment baaed upon the acceptability of the flying characteristics at 
high angles of attack. The end point descriptors selected by the pilots 
were, in their jargon, "Sierra Hotel," meaning they would like to fly 
the aircraft at high AOA and "Delta Sierra," meaning unacceptable flying 
characteristics. The resulting rating scale is shown in Figure 73. 
Both test pilots were happy with this final form and considered that it 
covered all key factors, minimized the descriptor conflicts, and did not 
constrain the pilot's quantitative evaluation. 

The rating scale of Figure 73 was employed in the final set of 
tracking tasks. Resulting characteristic assessments for the six 
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Complexity Required Timing 



Figure 73- Modified Departure Rating Scale (Mod II) 






vehicle configurations are summarized In Figure 74. The numerical 
ratings shown for each pilot are the averages of separate ratings given 
for each of the three tracking tasks (T^ through T3). Considering the 
small statistical run sampling and the fact that all configurations were 
intended to have basic handling deficiencies, the results are considered 
to be very gocd. It Is to be expected that all ratings would tend 
toward the poor (high numerical value) end of thj scales and they gener¬ 
ally are. Most importantly, they also reflect the specific dynamic and 
handling variations intended to be exhibited between configurations and, 
in most Instances, the ratings of RC and JF demonstrate agreement in 
variational differences if not in absolute levels. The ratings of pilot 
RH tend to vacillate between those o t RC and JF and also (except for 
recovery control timing) tend to ague with the poorer rating. ThiB 
downrating probably reflects pilot RH's relative inexperience in depar¬ 
ture/ ’-pin/recovery and the violent PSGs which high-performance fighters 
often exhibit. Thus the ratings of RC and JF are considered the more 
significant. 

The separation between RC and JF ratings on warning and motion 
severity is consistent with differences in piloting techniques. Gener¬ 
ally JF is about one rating point harsher than KC because of a more 
rapid pull into departure and spin. This will be discussed in more 
detail later. 

Recovery control timing assessment shows very good agreement between 
the pilots except for Configuration C2. This primarily reflects timing 
of control release . The differential is about one rating point or less 
except Configuration C 2 where RC had difficulty with the two different 
spin modes the aircraft could have. Timing of control release deter¬ 
mined whether an oscillatory (recoverable) or flat (unrecoverable) spin 
was obtained. 

Time to recover and overall rating show almost identical trends and 
thus reflect predominant concern for loss of altitude in any departure/ 
spin situation. Configuration had an oscillatory, recoverable 
spin. Configuration D did not spin but could end up with appreciable 
altitude loss if the pilot kept fighting the "wallowing" departure mode 
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Control Complexity 








(described in Subsection II.B). All other configurations had flat spin 
modes that were difficult if not impossible to stop. The significant 
disapproval of the Configuration D pitch-up is also evident for pilot RC 
on both of these ratings. 

The rating scale for recovery control complexity may need additional 
refinement. Pilot JF rated all configurations at 1 because he merely 
released controls to initiate recovery and, by definition, this is as 
simple as things can be. Pilot RC, with less severe departures, was 
able to discern some influence of the Btick and rudder in effecting 
recovery and thio was also reflected in his being more critical of 
recovery control timing. 

However, aspects pertaining to recovery (control usage and recovery 
time) for configurations other than D are not necessarily significant 
because the simulation cannot be considered valid for spin and spin 
recovery. The motions were generally so violent and the out-the- 

windscreen sky/earth display so featureless that the pilots had to 
resort to the all-attitude ball for recovery. This ia not completely 
unnatural for the military pilots, but did make the task more diffi¬ 
cult. Additionally, the aerodynamic coefficients employed in the 
simulation do not represent steady spin conditions. Probably the most 
that can be attributed to recovery evaluations is that they exercised 
the rating scales. 

In summary, the departure rating scales of Figure 73 were accepted 
and auppurted by the pilotB and, baaed upon this small statistical 

sample, the numerical ratings given reflect characteristics "designed 
into” the vehicle configurations. The spread in ratings between the Air 
Force test pilots (generally about O.S on a 5-point scale) is consistent 
with spreads normally expected in Cooper-Harper ratings (about 1.0 on a 
10-polnt scale). Greater spreads generally reflect sensitivity to 
piloting technique; that sensitivity itself may be an attribute of the 
scales. 

Since vehicle attributes rather than a task performance level 

are being rated, use of the rating scales must be accompanied by a 
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qualitative assessment/description. The qualitative information should 
include, as a minimum, the following: 


• Warning 

- Type 
Clarity 

- Margin 

• Departure 

« Resistance (susceptibility) 

- Type 

- Severity 

- Ability of pilot to delay or prevent 

• Control action taken 

e Demands on the pilot 

• Post-Departure Motion 

- Type of aircraft motion 

- Severity 

• Recovery 

- Rapidity 

- Recovery controls 
• Demands 

« Ability to recognise 

a Ability to perform necessary control action 


B. CONFIGURATION ASSESSMENTS 


As noted previously, the goals were to validate the influence of the 
key aerodynamic coefficients in determining departure characteristics} 
evaluate the influence of varied maneuver-limiting factors on high AOA 
maneuvering control; and identify potential flying quality criteria 
in terms of departure resistance, warning, and severity for possible 
application to the flying quality specification. In this subsection, 
predicted vs. observed characteristics aie reviewed and demonstrate 
that the above goals have been achieved, although there were seme sur¬ 
prises in the process. The unaugmented airframe configurations will be 
assessed firot, and then the augmented cases. 






It Is necessary to rely heavily on qualitative statements contained 
in the debriefing questionnaire and commentary taped on-line as the 
pilots talked themselves through the runs. The questionnaires were 
filled in only during the initial progression through the test matrix. 
At that time the pilots were "feeling out” the various configurations 
and produced voluminous commentary on each, covering the primary obser¬ 
vations plus variations due Co different control applications, tech¬ 
niques, rates of onset, tracking maneuvers, etc. The assessments 
reported in the following paragraphs have been distilled to the simplest 
possible factors consistent with the goals of the simulation. 

Only assessments of the fighter test pilots are reported because RC 
and JF are most experienced in stall/departure/spin testing and were 
more observant of differences between configurations. It became 
apparent early in the data analysis that the two test pilots were 
employing almost opposite approaches to stall/departure. RC employed a 
cautious, slow increase in AOA and generally was able to detect subtle 
changes in vehicle stability or response characteristics. JF used an 
aggressive, rapid AOA Increase which did not give time to detect such 
warnings. As a consequence, he generally pulled to a higher AOA with 
leas speed bleedoff and obtained more severe post-stall gyrations 
(PSG). Thus, differences in pilot closed-loop control techniques pro¬ 
duced a confounding Influence which must be recognized before delving 
into other details! departure warning and severity are a function of 
pull-up rate and lateral-directional controls application immediately 
preceding departure. 

1. Unangamated Airframe 

a. Configuration Aj 


Predicted versus observed pre- and post-departure characteristics 
are shown in Table 10. The warning Pilot RC observed with increasing 
AOA Is consistent with the graceful (gradual) degradation suggested 
by the open-loop dutch roll and roll numerator root migrations of 
Figure 55. It appears the negative C m ^ contribution was observed as a 
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TABLE 10. CONFIGURATION A, MANEUVER-LIMITING FACTORS 


PREDICTED 

OBSERVED | 

RC 

JF 

Warning: 

Wing rock 

(16 < o < 22 deg) 

Wing rock low ampli¬ 
tude, high frequency 

None 

Roll reversal 
(a > 18 deg) 

Roll reversal 

(a > 20 deg); g-break 

(20-22 deg a) 

None; g-break 


Wing rock large ampli¬ 
tude, low frequency 


Rating: 

2.8 

4.8 

Departure: 

Noee slice 
(o > 22 deg) 

Nose slice followed by 
roll 

Nose slice and roll 

Post stall gyrations 

Large yaw, pitch, roll 
oscillation and spin 

Violent oscillations 
and spin 

Severity Rating: 

2.2 

4.6 

Departure/Spin 
Susceptibility 
(Weissman criterion): 

a < 21 deg: mild or 
no departure 

Release a - 20 deg; 
yaw initially slow, 
builds rapidly 


21 < a < 24 deg: 
moderate rolling 
departure 


Release a ■ 23 deg; 
abrupt nose slice, 
builds exponentially 

a > 24 deg: strong 
rolling departure, 
high spin suscepti¬ 
bility 



Rating: 

Resistant 

Extremely Susceptible 

Comments: 

Maneuvers tend to mask 
warning; not masked if 
pay attention 

Moat violent post 
stall gyration oscil¬ 
lations of any con¬ 
figuration 
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"g-break." The nose slice departure Is as predicted, but the onset is 
noticed prior to C n ^ becoming negative. Thus onset might be influ¬ 
enced by pilot control, and thuB associated with negative 1/T^. Pilot 
RC considered the configuration departure resistant (R) because of ample 
warning and his tendency to initiate recovery at a * 20 deg where the 
Welssman criterion predicts mild or no departure. Pilot RC therefore 
gave fairly good high AOA ratings. 

Pilot JF eaw no warning except the "g-break." He obtained an abrupt 
noBe slice and violent post-stall gyration (PSG). The tendency to flat 
spin may be due to the large negative C np and small 0^^ at his recovery 
initiation point and the further destabilizing influence of sideslip on 
the denominator and numerator roots, all of which portend high body-axis 
yaw rate and little roll. Inertia effects would tend to raise the nose 

and align the yuw rotational axis with the velocity vector-a flat 

spin. 

This configuration also exhibited violent oscillatory spin charac¬ 
teristics which could be reduced if a large nose-down pitch rate devel¬ 
oped as the nose was slicing. This could be produced via forward stick 
coinciding with the "g-break." The resulting rq inertia cross-coupling 
produced a large roll acceleration and quite wild gyrations (s.g., 
Figure 75). Pilot JF noted this configuration to have the most violent 
PSG characteristics of any configuration and rated it extremely suscep¬ 
tible (ES) to departure. 

While the Weiaeman criterion correctly predicted the spin suscep¬ 
tibility, the initial departure was in yaw rather than roil. The 
maneuver-limiting factors are the strong nose slice and Bpin character¬ 
istics. The causal factors are combined open- and closed-loop (roll) 
instabilities accompanied by aignifeant negative C^. 

b. Configuration B 

As predicted, this configuration is similar to Configuration A, but 
with less departure warning (see Table 11). The low-frequency yaw/roll 
oscillation that Pilot RC observed above 21 deg AOA is consistent with 
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Figure 75. Configuration Aj Departure and PSG from Wings Level 
Tracking Taak; Pilot RC 
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TABLE 11. CONFIGURATION B MANEUVER-LIMITING FACTORS 


PREDICTED 

OBSERVED | 

RC 

JF 

Warning: 

Roll reversal 
(a > 18 deg) 

Resists roll command; 
g-break 

Low frequency yaw/roll 
21 deg a 

No warning 

Rating; 

4.2 

4.8 

Departure! 

Nose allca 
(a > 22 deg) 

Nose slice 

Nose slice 

Poat atall gyratlona 

Violent roll reversal 
or spin 

Severe and violent 
oscillation 

Severity Rating! 

2.7 

4.5 

Departure/Spin 
Susceptibility 
(Weissuan Criterion) 

a < 21 deg: mild or 
no departure 

Release a * 21 dag; 
yaw followed by roll 
and pitch down 


21 < a < 24 deg: 
moderate rolling 
departure 


Release a a 24 deg 

a > 24 deg! strong 
noae slice, high spin 
susceptibility 


Nose drift followed by 
roll; prone to flat 
spin 

Rating! 

Susceptible 

Extremely Susceptible 

Comments; 

Cannot discern onset 
of departure; insidi¬ 
ous 

Low acceleration In 
noao slice 
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open-loop dutch roll dynamics (see Figure 57). The lower frequency 
resulted In less acceleration Into the nose slice. Coupled with roll 
reversal, this apparently masks the departure onaet. Although he Ini¬ 
tiated recovery at about the same AOA and thought the departure severity 
(pilot rating ■ 2.7) Is about the same as Configuration Aj, he consid¬ 
ered Configuration B to be departure susceptible (S) because of the poor 
warning (PR ■ 4.2). Pilot JF could see little difference between A} and 
B except the increased tendency to flat spin due to the combined high 
aerodynamic roll rate damping and strong directional divergence. He 
again rated it extremely susceptible (ES). 

This configuration exhibited all anticipated flying characteris¬ 
tics. The Welssman criterion correctly identified departure and spin 
severity but not type. Again, the maneuver-limiting factors are the 
severe nose slice and flat spin tendencies. The causal factors are the 
open- and closed-loop directional Instabilities aggravated by high 
negative C£ p . 

c. Configuration 


Configuration Cj (Table 12) with positive dynamic stability (due to 
large negative C*^) had a much greater tendency to develop roll oscilla¬ 
tions. Both pilots noted that any lateral control input excited wing 
rock and produced a tendency to lateral PIO. Since this sensitivity 
persisted over a wide AOA range, it was not considered a warning but 
actually a masking of departure onset. Both pilots considered the 
warning to be poor (PR > 4). 

Thle configuration Is stable open-loop and driven unstable and Into 
departure by closing the attitude loops. Departure characteristics are 
less consistent (predictable) because the nature of departure depends 
upon pilot control activity. It exhibited at least two and possibly 
thrse different types of departure onset, which further aggravated the 
poor warning. Thu lateral phugoid mode could contribute the additional 
departure mode (see Figure 59). Ons pilot detected mors tendency to yaw 
and the other more tendency to roll. However, they generally were in 
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TABLE 12. CONFIGURATION C, MANEUVER-LIMITING FACTORS 



OBSERVED | 


RC 

JF 

Warning: 

Wing rock 

(16 < a < 26 dag) 

Wing rock: any 
disturbance 

Wing rock: tendency 
to low frequency 
pilot-induced 
oscillation 

Roll reveraal 
(a > 13 deg) 

Low frequency, high 
amplitude oscillation 
masks departure 


Rating: 

4.3 

4.0 

Departure: 

Roll departure 

Post-stall gyrations 

Nose slice 

Roll 

Two distinct types: 

a) roll, little yaw; 

b) yaw, followed by 
roll 

Severity Rating: 

1.6 

3.7 

Departure/Spin 



Susceptibility: 

(Welssman Criterion) 

a < 21 deg: no depar¬ 
ture 



a > 21 deg: mild 
rolling departure, low 
spin susceptibility 

Release a * 32.5 deg, 
nose slice and roll 

Release a * 25.5 deg, 
milder departure, 
oscillatory spin 

Rating: 

Susceptible 

Susceptible 

Comments: 

Takes time to see yaw 
is not going to stop 

Not prone to flat spin 









closer agreement for this configuration concerning warning, nature, and 
severity of and susceptibility to departure, probably due to the posi¬ 
tive dynamic stability which prolonged the time period the airplane 
could be maintained under semi-control. Thus, both pilots had suffi¬ 
cient time to observe its high AOA characteristics fully before closed- 
loop instability predominated. 

Both pilots considered the departure and PSG severity to be somewhat 
lass than for either A^ or B. The configuration was spin-prone, but 
spins were generally highly ocillatory with some possibility of recov¬ 
ery. It was not prone to flat Bpins. Departure and spin characteris¬ 
tics and susceptibility are in agreement with the Welssman criterion. 

The initial or basic maneuver-limiting factor is the strong wing 
rock tendency which extends from relatively low AOA through departure 
and is aggravated by closed-loop instability (roll PIO). Departure is a 
secondary factor. Warning may also be adversely influenced by the low- 
frequency, large-amplitude, lateral phugoid oscillation. 

d. Configuration D 

This configuration (Table 13) had less directional instability but 
did have sideslip-induced pitch-up. The expected wing-rock tendency 
from a lightly damped dutch roll was not observed by either pilot. Both 
noted departure onset as a slow, persistent yaw (nose slice) accompanied 
by mild roll reversal and pitch-up. The onset was sufficiently mild 
thet departure was semi-controllable throughout the achievable AOA 
range. A poor departure warning rating (PR > 3.5) resulted. 

Open- and closed-loop instability parameters and u)2 or C n ^ 
and LCDP) are small for this configuration and lead to sloppy control 
and low-frequency wallowing which masks daparture. At the same time, 
tha wallowing does not generate sufficiently rapid motion to excite 
inertia cross-coupling and PSG. All pilots tended to continue fighting 
to maintain control well past full stall, incurring excessive altitude 
loss. However, If controls were released st any tims the aircraft would 
Immediately go into a nosa-low spiral sod recover by itself. 
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TABLE 13. CONFIGURATION D MANEUVER-LIMITING FACTORS 


PREDICTED 


RC 


Warning! 

Wing rock 

(16 < a < 26 deg) 


OBSERVED 


JF 


Roll reveraal 
(o > 20 deg) 


Rating : 


Departure : 


Yaw, roll reversal, 
pitch up 

Sloppy control, large 
6 and $ overahootg 
mask departure 


Large definite yaw 
followed by reveraal 
and pitch up 


4.0 


3.5 


Pitch up 


Note alice/pitch up 


Noae alice/pitch up 


Poet atall gyration 

Departure/Sp in 
Susceptibility 
(Walasman Criterion): 


None, never apln 


Mild, no spin 


All a: no departure 


Release a * 27.6 deg, 
slow noae slice, pitch 
up, some roll 


Resistant 


Comments: 


Do not like pitch up 
at departure; cannot 
get noae down 


Release a * 27.6 deg; 
alow conalatent nose 
alice, roll, departure 
nose high 

Resistant 


Always recoverable If 
release controls 








Departure warning was considered poor (e.g., 4.0) by RC because of 
the wallowing, but at the sane tine the configuration was considered 
departure resistant (R). Pilot JF assessed It as resistant primarily 
because it would not spin. Pilot RC's principal objection to this 
configuration was the positive C^. Even though it eliminated the 
severe pitch-down in nose slice and the resulting inertia-coupled PSG 
exhibited by other configurations, he did not like the nose rising, no 
matter how mildly, when in or near departure conditions. 

The maneuver-limiting factor for this configuration is the wallow¬ 
ing - a general insensitivity to or sloppiness of response to pilot 

control inputs prior to the departure AOA. The key causal factors are 
considered to be the relatively small (due to small C n0 and Cj^) 
which produces very low DC gain in the roll frequency (Bode) domain and 
positive C m p which apparently contributed to overcontrol in pitch. 

e. Summary 

The intended high AOA departure characteristics were obtained in the 
four unaugmented airframe configurations. The results demonstrated the 
iniivences that changes in the six key aerodynamic coefficients ctn have 
on departure warning, severity, and susceptibility. Tht .susceptibility 
of Configuration Cj^ to departure was somewhat revealing! the open-loop 
parameter C n ^ , by itaelf, is not a sufficient indicator of departure 
resistance. However, when combined with the closed-loop parameter LCDP 
(or «J), the two are quite accurate indicators of departure suscepti¬ 
bility and aeverity as advanced by Weissman. 

Departure severity also appears Influenced by the static cross¬ 
coupling derivatives £&, Nq, and Mg. Severe departure and PSG were 
obtained for the two configurations in which combinations of these 
derivatives shifted lateral-directional denominator and numerator roots 
into the right half-plane of the root locus plot with M 0 (Figures 55 
and 57). Less severe departure was obtained for combinations which 
primarily shifted numerator roots into the right half-plane (Figure 59), 
and the mildest departure was obtained for combinations which produced 
very little shift of lateral denominator and numerator roots with 
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sideslip (Figure 61). The literal approximate factors developed In 
Section III Indicate that the small shifting Is due to the combined 
Influence of small N£ and positive Mg. However, the nose rise asso¬ 
ciated with large positive Mg was disconcerting to one pilot. 

The simulation results also demonstrated that wing rock or lateral 
PIO Is not a satisfactory warning of Impending departure if '■.here is 
sufficient elevator power to pull through the warning region rapidly or 
if the AOA range In which it occurs is too wide. 

2. Augmented Airframe 

a. Configuration A2 


The flight control augmentation was configured to provide a high 
level of roll rate damping, and to improve roll maneuvering response via 
roll rate command (p c ) and reduced adverse aileron yaw (SRI). This 
augmentation was also expected to eliminate all cues of impending depar¬ 
ture. The departure was anticipated to be a strong nose slice due to 
increasingly negative C n0 as AOA increased and, once sideslip started to 
build, it was expected that the characteristics would be very similar to 
those of Configuration 

Results summarized in Table 14 for Pilot JF are very much as antici¬ 
pated. However, Pilot RC detected a pre-nose-sllce resistance to roll 
commands (at AOA above approximately 21 deg) when attempting to track 
the maneuvering target. This resistance was observed only during the 
initial tracking sets and resulted in the pilot's neutralizing all con- 
trola before the directional instability became severe (rating 1.8). 
The aecond set of tracking runs were done mure hurriedly, the warning 
waa not obtained (rating 4.3), and more severe departures ensued (rating 
3.7). The susceptibility assessment as rssistant (R) by Pilot RC was 
given only on the initial set of runs. 

The maneuver-limiting factor for this configuration is the nose 
■lice. But. since there la no warning, an AOA limit would have to be 
set to provide a safety margin against inadvartant departure. The key 






146 







PREDICTED 


OBSERVED 



Warningt 

None 

Resists above 

21 deg a 

Claar warning with 
target 

None 

Rating; 

2.8 

<4.3 on reruns) 

4.8 

Departure; 

Nose slice 

Nose slice 

Nose Blice 

Roll and severe PSG 

Severity Rating 

1.8 

(3.7 on reruns) 

4.5 

Departure/Spin 
Susceptibility 

o < 23 deg; no depar¬ 
ture 



a > 23 deg: nose 
slice and revert to 
Configuration Aj 
characteristics 

Release a * 25 deg; 
Bmooth, slow yaw 
followed by roll 

Release a 4 32 deg; 
exponential yaw diver¬ 
gence tendency to flat 
spin 

Rating; 

Resistant 

Extremely Susceptible 

Comments: 

Depart without pilot 
aware; only cue is 
resistance to $ c 
against target 

Easier to fly 

Rate of yaw dependent 
upon 4 or g 









causal factors are the strong static directional divergence and the 
elimination of departure warning by the augmentation system. 

Time traces of a Configuration A 2 departure from a wings-level pull- 
up are shown in Figure 76. The lateral Btick trace reflects about the 
same activity prior to departure as ir. Figure 75 but far leas sideslip 
and roll oscillation. After pilot neutralization of controls, the aug¬ 
mentation loops continue to deflect the lateral and directional control 
surfaces and apparently aid recovery (again compare with Figure 75). 
Interestingly, the actual departure and PSG motions of Figure 76 are 
remarkably similar to Reference 33 flight traces for the F-4E with 
augmentation systems on (Figure 77). Prior to departure the traces of 
the two figures are not comparable because the stick to rudder intercon¬ 
nect in the simulation reduced sldoslip excitation and provided better 
roll attitude control. 


b. Configuration C 2 


Based on the positive C ns and improved w? (or LCDP), this aug- 
Pdyn “ 

mented configuration was expected to have the best high AOA flight 
characteristics and little or no departure tendency. It did allow 
both pilots to consistently maintain control to higher AOA than did any 
other configuration, though abrupt nose slice departures ware common 
(Table 15). In this case the roll damper countered the dynamic 
stability contribution and allowed the largo negative otatic insta¬ 
bility to predominate. The aircraft responses by themselves provided no 
warning of impending departure, although Pilot RC was able to extract 
some low-frequency PIO cues with respect to the target aircraft. He 
rated the configuration doparture resistant (R) because of this cue and 
the mild (to him) nose olice motion. Without ths target he rated the 
configuration departure susceptible (S), as did Pilot JF. 


This configuration was also spin-prone and, similar to Cp exhibited 
two distinct spin modes as noted in Table 15. It appeared to be more 
prone to spin flat than was the unaugmented configuration (C t ), because 
the roll rate CAS automatically opposed any uncoromandad roll. 
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Figure 76. Configuration A2 Departure and FSG from Wings Level 
Tracking Task; Pilot RC 
























Figure 77. Typical F—4E Roll Departure/Spin Entry 














TABLE 15. CONFIGURATION C, MANEUVER-LIMITING FACTORS 



OBSERVED | 

RC 

JF 

Warning: 



None 

Low frequency pilot- 
induced oscillation; 
g-break 

None, wry steady 

Rating: 

3.2 

4.8 

Departure: 



None 

Nose slice followed by 
roll and pitch down 

Abrupt nose slice 

Severity Rating: 

1.5 

3.8 

Departure/Spin 

Susceptibility: 

Resistant 

Release a * 31 deg; 
slow nose slice; 
almost always spin 

Release o * 30 deg; 
abrupt nose slice; two 
spin modes: 

a) oscillatory 
(recoverable) 

b) flat (nonrecover- 
able) 

Rating: 

Resistant 

Susceptible 

Comments: 

Resistant because tar¬ 
get gives cue; without 
target is Susceptible 

Better flying quali¬ 
ties prior to depar¬ 
ture 









The maneuver-limiting factor for this configuration is the nose 

slice departure. The positive C ni9 and throughout the achievable 
Pdyn v 

AOA range did allow it to be controlled to considerably higher AOA than 
the other configurations. In this respect it was less susceptible to 
departure. However, at the higher AOA directional Instability even¬ 
tually required full opposing rudder deflection, which was immediately 
followed by departure. 

c. Overall Assessment 


The augmented flight control system designed to enhance character¬ 
istics in normal flight improved pre-departure control precision, 
reduced departure warning, and increased the tendency to flat spin. The 
roll rate CAS and the SRI together improved lateral-directional control, 
which allowed higher AOA to be reached before sideslip became signifi¬ 
cant. The major portion of thiB benefit derived from the SRI. For 
Configuration C2 the roll rate damper suppressed the <£g$ contribution 
to dynamic stability and thus increased the nose-slice departure and 
flat spin tendencies. Thus, this augmentation function actually 
degraded high AOA characteristics. 

C. SUMMARY 


The intended high AOA departure characteristics were obtained in the 
four bare-alrframe and two augmented configurations. The results demon¬ 
strate the strong Influence changes in the six key aerodynamic coeffi¬ 
cients (£a> «£p» Na» and Mg) can have on departure warning, 

severity, and susceptibility. The simulation plus literal approximate 
factors of Section III show: 

• The most violent departure characteristics and 
flat spin tendem. Lea to be associated with large 
negative N^, Ng, and Mg. 

• The least violent departure characteristics and 
greatest spin resistance to be associated with 
small to moderate negative and Ng and positive 
Mg- 
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• The tendency tc Inconsistent departure/spin char¬ 
acteristics (possibly several modes) to he asso¬ 
ciated with positive Cngjy n brought about by 
large negative C^.. and continued large adverse 
yaw <e»|, LCDP). P 

However, the observed high AOA stall, departure, and spin dynamic char¬ 
acteristics also depended significantly on pilot control technique. 
This is consistent with actual flight experience in which It haB been 
observed, for example, that spin characteristics and recovery techniques 
can be quite different, depending upon the maneuver from which the spin 
occurred. 

The final high AOA departure/spin rating scale of Figure 73 was 
accepted and supported by the pilots. Results of this first application 
showed consistent reflection of both the airframe dynamic characteris¬ 
tics and any sensitivity to difference in piloting technique. However, 
the ratings must be accompanied by detailed qualitative description of 
the characteristics observed, i.e., what is being rated. 

In this vein, it might be noted that there are some differences 
between ratings for the initial familiarization and tracking runs 
reported in Subsection B and for the final tracking runs reported in 
Subsection A. These are due to several factors. First, Pilot RC was 
working with the original open-ended scale (Figure 72) during the ini¬ 
tial runs and the cloaed-ended scale (Figure 73) during the reruns; 
Pilot JF used the close-ended scale for both sets of runs. Second, 
during the first sequence of runs Pilot RC was unhurried in investi¬ 
gating the characteristics of each configuration. However, for the 
repeat runs he was under some pressure to complete the series of maneu¬ 
vers and configurations within time constraints. Thus, the rate at 
which he approached stall/departure was more rapid and this influenced 
his assessment of warning and motion severity, as demonstrated in 
Table 16. The greatest differences are for the augmented flight control 
configurations, where RC saw considerably less warning and more severe 
departure during the reruns. Results for JF, shown for comparison, 
reflect quite consistent ratings between the sets, since he always used 
a rapid pullup technique. 
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TABLE 16 


COMPARISON OF PILOT RATINGS 
INITIAL VS. RERUNS 



CONFIGU¬ 

RATION 


INITIAL RERUN 


MOTION SEVERITY 


INITIAL RERUN 


MOTION SEVERITY 


INITIAL RERUN INITIAL RERUN 










SECTIOI VII 


DEPARTURE SOSCEPTIBLITT (XITERIOB 

Pitch control techniques and aggrasalveness had obvious Influences 
on the departure warning and susceptibility assessments shown In the 
previous section. When asked about specific cues and personal criteria 
which would prompt the pilot to acknowledge departure and therefore 
to Initiate recovery, vague answers were obtained, e.g., "persistent 
or continued uncommanded roll or yaw" or "sudden, rapid uncommanded 
motion — somewhere around 20-30 deg/sec.” However, detailed scrutiny 
of the time traces failed to reveal correlation between release of back 
stick and various motion quantities (p, r, \p, 0, etc*). These quan¬ 
tities were always changing rapidly and timing became a big factor. 
Deviations of fractions of a second In Initiating recovery (or reeding 
the tracao) produced a large variation in the motion amplitudes (e.g., 
see Figures 57 and 58). The greatest correlation was found In AOA 
Itself — possibly because Its rate of change was the slowest of all and 
therefore could be read and correlated with the greatest accuracy with 
what the pilots saw. 

Since the piloting task required closed-loop attitude control,, and 
there was considerable evidence of LCDP or being a key parameter, 
attention was turned to the possibility of a closed-loop criterion. A 
simple analysis of theoretical roll loop stability was found to agree 
quite well with frequencies, damping ratios, and divargences obtained 
from the time traces. Furthermore, the results corroborated the differ¬ 
ences In pilot techniques and observed warning reported in the previous 
section. 

This section summarise# the relationship obtained between the 
closed-loop analysis and the piloted ulraulatlon time traces and presents 
a closed-loop criterion which is consistent with and explains the seem¬ 
ingly diverse departure susceptibility assessments given by the two Air 
Force test pilots. 
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A. CLOUD-LOOP ANALYSIS 

The lateral stick time trace of Figure 57 (Configuration A^) le 
typical of most runs for Pilot RC. The phasing between 4 excursions and 
stick deflections Indicates little, if any, lead being generated during 
the PIO at a - 16 deg and above. Therefore, the siaple block diagram of 
Figure 78 was chosen to model the unaugmented airframe, control system, 
and pilot. A Bode-root locus closed-loop survey plot for Configuration 
Aj at a ■ 16 deg and 0-0 deg is presented in Figure 79. Closed-loop 
gains and frequencies obtained from strip chart recordings of simulation 
runs for Pilot RC are shown on the Bode amplitude plot by the eeveral 
horisontal lines (K^) and tic marks (-+■)• [For example, at 16 deg AOA 
the traces of Figure 57 show a pilot gain * - 4.5 (13 dB) and fre¬ 
quency w * 1.8 rad/sec.] The frequency of the undamped oscillation 
obtained from the strip chart is in remarkable agreement with the fre¬ 
quency for 180 deg phase of the analytical predictions, validating the 
assumptions that this pilot is adopting neither lead nor lag in attempt¬ 
ing to control roll attitude. Rather, he la merely reacting with a 
normal time delay with stick proportional and opposite to bank angle 
excursions. 



















The measured gains range from about 10 to 20 dB. This uriatlon. 
Identified on the root locus plot of Figure 79, Indicates but small 
Influence on dutch roll divergence. The 13 dB gain point Indicates a 
mild divergence which is also in agreement with the traces of Figure 57 
at 16 deg AOA. 

As AOA is Increased in Figure 57, the lateral stick amplitude 
remains fairly constant but the <j> excursions increase. Thus, the loop 
gain (K^) decreases slightly, the oscillation period increases to about 
6 seconds (w * l rad/sec), and a yaw rate divergence develops. Again, a 
closed-loop survey plot for a ■ 23 deg (Figure 80) shows excellent 

agreement, l.e., uij ■ 1 rad/sec and a firot'-order spiral divergence 
of T; * 2 sec which results from the spiral mode being driven toward the 
RHP zero 1/T^. 

Similar trends were obvious from Pilot RC run traces for moat con¬ 
figurations. He appeared to be following the neutral stability boundary 
of the oscillatory mode and initiating recovery when the first-order 

divergence mode exceeded some as yet unidentified divergence rate. 

Therefore, using the same pilot model, the maximum boundary for 

stable modes was calculated and plotted as shown in Figure 81. This 

shown two boundaries. The upper la the gain not to be exceeded at B « 0 
and the lower is the reduced limit gain for 3 » 5.S deg. To the right 
(higher AOA) of either boundary a lateral or directional instability 
exists and if AOA is increased must be decreased to maintain sta¬ 

bility. 

The roll loop gains employed by both pllotB were determined from the 
time traces and located on similar stability boundary plots drawn for 
each vehicle configuration. Figure 82 is the plot for Configuration A^ 
again, but now showing gain bands for each pilot. The band for RC shows 
that he tended to keep his gain within 10 dB above the maximum stable 
boundary and to decrease his gain as AOA was increased. On the other 

hand, JF started with a lower gain and tended to keep it constant during 

the pull-up. 
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Figure 82. Roll Loop Closure Gains Employed with 
Configuration 







Figure 83 Is the gain plot fot Configuration Aj • Again JF used a 
lower, constant gain while RC tended to follow the stability boundary 
until the sudden drop. Figure 84 is the plot for Configuration C^, 
showing results similar to thoiie for Configuration Aj. This configura¬ 
tion Is somewhat unusual In that a second stable region exists for 2^0 
at very high AOA. Figure 85 for Configuration C 2 is similar to Configu¬ 
ration A 2 except there is no dropoff in stability boundary at high 
AOA. The results for Configuration D (Figure 66) are consistent: Pilot 
RC followed the stability boundary and hence should detect departure 
onset; Pilot JF used a lower, constant gain which suddenly crosses the 
boundaries as AOA is Increased. He therefore obtained little warning. 
The plot for Configuration B is presented in Figure 87. This vehicle is 
so well damped In roll that it was difficult to obtain the necessary 
data from time traces except for one run by RC. 

At the bottom of Figures 82-87 negative values of the open-loop roll 
numerator root 1/T^ are plotted versus a (at 8 - 0). Figure 85 does 
not show this extra plot because wj- 1 b positive to o > 35 deg for this 
configuration. Also identified on each plot is the average AOA at which 
recovery was initiated by Pilots RC ( ^symbol) and JF (<?). It will be 
noted that: 


• In most instances a decrease In the boundary 
is related to l/T^ becoming negative (RHP aero). 

• In all but one case recovery is initiated by both 
pilots after exceeding the K, boundary and exper¬ 
iencing negative l/T^, (the* only exception ie 
Configuration C 2 for which the numerator roots 
always lie in the LHP and thers is no dropoff in 

y. 

• Pilot RC always employs higher roll loop gain, 
alwaya follows the K* boundary, and generally 
initlatee recovery at* lower AOA. However, for 
Configuratlona and C 2 he did explore higher 
AOAs before recovery. 

• Pilot JF tended to adopt a lower roll loop gain 
and kept it constant during his more rapid 
(aggressive) pull-up. He generally saw little or 
no warning and penetrated further into instabili¬ 
ty regions before initiating recovery. 
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Figure 85. Roll Loop Closure Gains Employed with 
Configuration C2 



Figure 86. Roll Loop Closure Cains Employed with 
Configuration D 
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Figure 87. Roll Loop Closure Gain Employed 
with Configuration B 


The foregoing observations appear to form the basis for each pilot's 
definition of departure susceptibility (cf. Section VI) s 

RC: Departure resistant (R) if warning is clear and 
consistent 

Departure susceptible (S) if warning Is unclear, 
masked, or Inconsistent 

JF: Departure resistant (K) 11: always recovered 

Departure susceptible (S) if sometimes recovered 

Departure extremely susceptible (ES) if never 
recovered 

Thus each pilot is viewing different aspects of the departure# Fortu¬ 
nately, this difference ierves to increase the information potential of 
the simulation program. This is probably typical of flight test, and 
indicative of why such sl:all/spin programs should be flown by more than 
one pilot. 
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B. POSSIBLE DCPAATUU SOSCBPTI1ILITT GKITRUA 

A comparison between depertute/epin susceptibility predicted by the 
Welsemen criterion end the assessments provided by the two pilots Is 
shown Is Table 17. Assessments substantially in agreement with predic¬ 
tion are shown in boxes. Obviously the aggressive pilot observed the 
worst possible characteristics of each unaugmented configuration as 
predicted by the criterion. The less aggressive pilot experienced 
something quite different. As noted previously, the criterion does not 
lend Itself to prediction of the frequency-dependent augmented airframe 
charcterlstlcs, but these configurations were expected to be less sus¬ 
ceptible to departure. This Influence was observed only by the less 
aggressive pilot. 


TABLE 17 

COMPARISON OF PREDICTED AND ACTUAL DEPARTURE 
SUSCEPTIBILITY ASSESSMENTS 


DEPARTURE SUSCEPTIBILITY 

CONFIGU¬ 

RATION 

PREDICTED 

RC 

JF 

A 1 

ES 

R 

H 

B 

ES 

S 

0 

Cl 

S 

0 

1S-ES 1 

D 

R 

0 

0 

a 2 

NA 

R 

ES 

c 2 

NA 

R 

S 


R - Resistant ES - Extremely susceptible 

S ■ Susceptible NA - Not applicable 
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One important difference between prediction and simulation was the 
nature of departure. Figure 63 predicts predominantly rolling departure 
with no indication of yaw departure. Our Configurations A, B, and D 
exhibited initial yaw excursions sometimes followed by roll. Configura¬ 
tion C exhibited two or even three different departure modes which were 
dependent upon control application at onset of departure. However, a 
rolling type motion did predominate. Thus, the Welssman criterion was 
not as useful in predicting the nature of departure observed by the 
pilot. 

The closed-loop analysis to Identify causal factors behind the 
widely differing pilot ratings for departure/spin susceptibility pro¬ 
duced a strong relationship between pilot ratings for both pilots and 
penetration of the roll numerator into the RHP. Figure 88 shows the 
value of the real part of the numerator root at the inetant the pilot 
decided he had departed, or was about to depart, and Initiated recovery. 
These values are plotted against the AOA at which recovery was started. 
(Note this is not the usual root locus ju> axis.) The points represent 
all six vehicle configurations as evaluated by both pilots. 

The division between departure-resistant (R) and susceptible (S) 
ratings is seen to lie at roughly -0.5 rad/sec. This corresponds to a 
time to double amplitude of approximately 1.4 sec. Zeros which lie to 
the left of this line apparently limit the first-order divergence to a 
rate slow enough for pilots to respond and recover. Zeros to the right 
of the line apparently allow divergence rates so fast that the pllotB 
cannot prevent departure. Again this Is consistent with the previously 
noted definition of departure whereby both pilots Indicated a threshold 
on rate of motion; however, the pilots were vague as to the value (e.g., 
"maybe 20 or 30 deg/sec"). One data point in Figure 88 violates the 
boundary. This is the augmented Configuration Cj. In this case the SRI 
eliminated adverse yaw and thus made the vehicle more departure resis¬ 
tant (as viewed by RC). However, it could be departed; and then the 
augmentation produced pro-apin control. Pilot JF rated apln suscepti¬ 
bility (and recovery) and apparently rated this configuration accord¬ 
ingly. 
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Figure 88. Departure Susceptibility Rating Versus 
Lateral Closed-Loop Divergence Potential 


The implication is that if the combined aerodynamics and flight 
control system design is such that 1/T^ never exceeds -0.5 throughout 
the achievable AOA range, the airplane will be departure-resistant. It 
should be noted that this criterion places no restriction on open-loop 

stability. For example, C n „ can be negative and, in fact, is nega- 
^dyn 

tlve for Configuration D (see Figure 62) which la rated dep.rture- 
realstant (see Table 13) by both pilots. 


Since u,* is the dimensional form of LCDP and since, in general, 
|1/T^|2 A lu^l* than one can relate the above 1/T^ boundary to an 
equivalent LCDP. For the flight conditions, inertias, etc., employed 
in this simulation, 1/T^ of -0.5 corresponds to LCDP of -0.001. This 
coincides with Welssaan's boundary between Regions A and B for positive 
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C no , see Figure 89: however, it is a little more conservative at 
^dyn 

negative C n . . Thus, the results of our simulation are compatible 
^dyn 

with and support Welssman's empirically derived LCDP boundary. The key 
difference in the criterion is that l/T^ Is not restricted to airframe 
lateral-directional static coefficients but can be applied to the com¬ 
pletely coupled 6 DOF airframe with a full complement of augmentation, 
stlck-to-rudder crossfeed, etc., throughout the aircraft development 
cycle. 


In addition, results of our simulation tend to indicate that another 
boundary might be appropriate in the upper left quadrant of Figure 89 to 
create two additional regions, E and F, as suggested in Reference 40. 
Region E would be classified as mild directional divergence and moderate 
spin tendency. Region F might be classified bb severe directional 
divergence and strong spin tendency. Note from Figure 63 that our air¬ 
craft Configuration D, which had mild departure characteristics but 
little or no spin tendency, would extend into Region E. 



Figure 89. Possible Modifications to Welssman Criterion 
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Additional support for the limit on 1/T^ can be drawn from the 
current MIL-F-8785C spiral divergence limit* The Level 3,* Category A 
and C Flight Phases requirement is that time to double amplitude be not 
lees than 4 sec* This limit came, in part, from the Reference 41 inves¬ 
tigation in which aircraft configurations having various spiral diver¬ 
gence races were rated by twelve pilots in a landing approach task* 
Results of that investigation, shown in Figure 90, indicate that diver¬ 
gences having a time to double amplitude of less than 2.4 sec were 
considered Intolerable. The time to double amplitude {T2) can be 
converted to an equivalent first-order time constant (T) by the rela¬ 
tionship 


T 2 * 0.74T 

A plot of pilot rating versus spiral-mode inverse time constant (1/T) is 
presented in Figure 91. The current -8785C Level 3 spiral limit is 
shown along with the proposed 1/T^ limit. 

The two limits are in good agreement when it is remembered that 
1/T a represents a pole (modal response) while 1/T^j represents a zero. 
Under closed-loop roll control the migration of an airframe pole to¬ 
ward 1/T^ would be dependent upon pilot gain; infinite loop gain would 
be required to obtain a closed-loop pole exactly at the zero location. 
A more realistic situation is shown in Figure 92, which is a Bode-aiggle 
plot for Configuration ^ it « ■ 21 deg, 0 ■ 0 deg. The dashed line 
terminating at 1/T^ - -0.6 is the locus of the closed-loop first-order 
divergence root for increasing gain, . Two typical gain lines for 
Pilot RC extracted from simulation traces at a - 21 ± 1 deg are shown. 
These closurea result in closed-loop roots at -0.2 and -0.33 rad/sec. 
Referring back to Figure 91, it may be observed that these straddle the 


Level 3: “Flying qualities such that the airplane can be con¬ 
trolled safely, but pilot workload is excessive or mission effectiveness 
is inadequate, or both. Category A Flight Phases can be terminated 
safely, end Category B end C Flight Phases can be completed." 








Rat,n Q o ? Average Pilot Rating 



Time to Double Amplitude of Aperodic Motion, T 8a (sec) 

re 90. Variation of Pilot Ratings with Time to Double Amplitude 
E Aperiodic Motion; Landing Approach Condition (from Ref. 40) 


-8785C 



Figure 91. Variation of Pilot Ratings with Equivalent 
Aperiodic Inverse Time Constants 


















1/T g boundary and are consistent with the Reference 40 division between 
tolerable and Intolerable handling. Again, our Pilot RC considered 
Configuration to be departure-resistant (R). 

The above discussion tends to show that the 1/T^ limit selected 
(-0.5) may not be conservative, since a very tight roll loop closure 
could result in closed-loop roots in the intolerable range of Figure 
91. Further investigation will be required to determine if the limit 
should be reduced. 


C. SUMMARY 

Analysis and piloted simulation have Bhown that: 

• The pilot's perception of departure suscepti¬ 
bility was found to be correlated with movement 
of one root of the roll numerator for lateral 
stick control into the RHP of the root locus, 
i.e., a nonmlnlraum-phaae zero. If pilots could 
fly the aircraft to such AOAs, they rated high 
AOA configurations which produced a zero, 1/Tx , 
more negative than -0.5 departure-susceptible. 

If this boundary was not or could not be ex¬ 
ceeded, the aircraft was considered departure- 
resistant. This rating is not a function of the 
sign or magnitude of the dynamic stability param- 

Cter Cn 0dyn* 

• A value of 1 /Ta. - -0.5 corresponds for the air¬ 
frame tested to 1 an effective LCDP of -0.001 and 
thus is consistent with and supports the empiri¬ 
cally derived LCDP departure boundary developed 
by Weiasman. However, results of the simulation 
were not in agreement with the types of departure 
predicted by Weiseman in that negative Cngdyn re ” 
glona produced yaw or nose-slice type departures. 

• A value of 1 /Ta. ■ -0.5 also is compatible with 
the current MIU^F-8785C Level 3 limit for spiral 
divergence in Category A and C Flight Phases. 

It is suggested that as the aircraft design/development cycle pro¬ 
gresses past the static wind tunnel phasr the Weiasman criterion for 
departure/spin susceptibility be replaced by a negative 1/T^ limit. 
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SECTION VIII 


OONOJOSIOHS AND RKGMfENDATIONS 

The goals of this research program were to! 

• Identify key design parameters that Umit high- 
angle-of-attack maneuverability for contemporary 
high-performance attack and fighter-type air¬ 
craft. 

• Postulate fundamental aerodynamic and control 
system design methodologies that will alleviate 
the limiting conditions. 

• Formulate handling qualities requirements for 
high-angle-of-attack maneuvering flight to be 
Incorporated in MII.-F-8785C, "Flying Qualities of 
Piloted Airplanes." 

It should be borne in mind that the results are based on variations on a 
single nonlinear aerodynamic model representing a region dominated by 
phenomena that ore highly configuration-dependent. Therefore, the gen¬ 
eralizations drawn are in need of further substantiation. 

Km KKT MANEUVER-LIMITING FACTORS 

The most crucial factors pertain to limiting the safe flight enve¬ 
lope, i.e., departure from controlled flight. By varying static aerody¬ 
namic coupling and cross-coupling coefficients, the nature of high AOA 
characteristics and susceptibility to departure were changed to reflect 
nose slice, rolling, and pitch-up departures. 

Nose slice (yaw) departures predominated for configurations having 
Cngdyn l&tge» negative. Analysis and prior simulation have shown that 
nose slice can be aggravated once sideslip startB by the static aerody¬ 
namic coefficients Cg^ and These aerodynamic coupling terms can 

cause a cross-coupled RHP zero to appear in the pitch numerator, ^6 sta ^* 
Any simultaneous attempt to control pitch attitude then produces or fur¬ 
ther aggravates directional divergence. The large negative 
figurations were prone to flat spins. 
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Rolling departures predominated when a large component made 

^ n 5dyn positive. These configurations also were prone to wing rock and 
lateral PIO. They exhibited at least two different PSG and spin modes, 
and were sensitive to maneuver and pilot activity at departure onset. 

Pitch due to sideslip, C m ^, of either positive or negative sign 
can adversely influence departure and handling characteristics. Posi¬ 
tive C m g can produce pitch-up which tends to mask a "g-break" stall 
warning cue and oppose pilot attempts to prevent stall/departure via 
lorward stick. Negative C m ^ tends h o augment "g-break" stall warning 
but reinforces pilot-commanded recovery pitching moment to the point 
that pitch-roll-yaw coupling may produce violent PSG. The presence of 
appreciable C,^ can be detected from flight traces of motion about a 
nominal zero Bideslip via a longitudinal oscillation at twice the dutch 
roll frequency. 

Unstable open-loop characteristic modes (poles) alone were not found 
to limit the safe flight envelope. If the pertinent control numerator 
zeros were favorably located (LHP), the pilot could prevent divergence. 
If the zeros were unfsvorably located (RHP) divergence could not be pre¬ 
vented; in this case the presence of open-loop instability increased 
departure severity. Thus, pole-zero locations had influence on depar¬ 
ture susceptibility and severity. In turn, pole-zero locations were 
shifted by static aerodynamic cross-coupling derivatives (X 3 , N^, Mg) 
when the magnitudes of these coefficients approached those of the common 
static stability derivatives («Cg, Ng, Mq): 

• c£ Ng, and Mg, in combination, influence the 
high-frequency modes of both the lateral and 
longitudinal denominators. 

• Nq and Mg, in combination, influence the zeros of 



• Nq and in combination, influence the zeros 
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B. ADDITIONAL FINDINGS AND OONOLOSIONS 


One major contribution has been identification and validation that 
pilot perception of lateral-directional departure susceptibility is 
related to one zero of the numerator becoming negative. Root 

magnitudes more negative than -0.5 rad/sec were consistently rated as 
departure-susceptible, while those less negative (or positive) are rated 
as departure-resistant. This criterion reflects a closed-loop diver¬ 
gence rate limit related to the pilot's threshold for uncommanded motion 
or ability to cope. As such it is a pilot-centered criterion which 
should be applicable for any flight situation, although it has been 
Identified in a low-Mach-number, fixed-base simulation. It is consis¬ 
tent with the empirically established airframe-alone departure/spin 
criterion boundaries of Weissman and extends applicability of that 
criterion to highly augmented airframe cases. It is also consistent 
with previous in-flight simulation of maximum controllable aperiodic 
divergence rates. Finally, it serves as both a design guide and a 
flying quality specification item. 

The “unstable" (RHP) zero generally results from negative (adverse) 
Cng and Cnfigg^. It therefore can be alleviated by aircraft configura¬ 
tion and flight control system oodifications which principally alter 
those coefficients. Ab noted above, the zero also can be highly influ¬ 
enced by the static aerodynamic cross-coupling coefficients and Cj,g. 

A second major contribution has been the development and partial 
validation of a flying quality rating scale applicable to the stall/ 
departure/recovery flight regime. It is believed that the scale has 
potential for defining the aircraft permissible flight envelope in much 
the same manner as the Cooper-Harper scale is now used in defining 
flying quality bounds. 

The simulation demonstrated the importance of loss-of-control 
warning in pilot assessment of high AOA flying qualities and departure 
susceptibility. A gradual and consistent degradation in stability and 
controllability was found to provide the best warning; however, the 
realization of such warning can be highly dependent upon pitch control 


•:,v 
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power. Oscillations such as wing rock or lateral PIO cannot be counted 
on for departure warning in air combat if the aircraft can be pulled 
through the "warning" region before the oscillation becomes apparent to 
the pilot. 

Low-frequency dutch roll or lateral phugoid oscillation was found to 
mask departure onset. Roll about the velocity vector at high AOA pro¬ 
duces significant "nose" yaw. The pilot has considerable difficulty 
differentiating between nose slice onset and low-frequency roll oscilla¬ 
tion about the velocity vector. High roll rate damping (natural or 
augmented) was also found to reduce high AOA departure warning and 
increase flat spin tendency. 

A flight control or augmentation system mechanization which improves 
and extends the aircraft controllability range by removing natural 
warning can actually be a detriment if the aircraft can still reach some 
departure AOA. To be of real benefit such a flight control system must 
also prevent departure. 

A stick-to-rudder interconnect (SRI) can be beneficial in reducing 
adverse "aileron" yaw, thereby providing more favorable LCDP or 1/T^ 
values. However, if static aerodynamic cross-coupling is strong, the 
analysis of Configuration A 2 has shown that even small sideslip can 
shift the numerator roots in a manner to negate the SRI contribution. 
These results also imply that a finite-authority S-llmiter may be of 
little benefit as a departure preventer if one or more RHP zeros occur 
in control loops the pilot is closing. That is, pilot control can drive 
the vehicle into a divergence which may exceed the (J-liraiter's capa¬ 
bility. 

Results of this study demonstrate that linear or quasi-linear 
frozen-point analysis is applicable to high AOA situations, but consid¬ 
erable caution must be exercised in the interpretation of results where 
aerodynamics are strong functions of angles of attack and sideslip. Any 
analysis pertaining to such regions must be based on a 6 DOF model with 
aerodynamic and kinematic cross-coupling terms, must Include non-zero- 
sides lip trim points, must view transfer-function numerator as well as 
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characteristic-equation parameters, and must be supported by nonlinear 
simulations. 

In particular, the results demonstrated the power of vector polygon 
and closed-loop root extraction analytic techniques for Identifying 
cause-effect relationships in highly coupled dynamic systems. The tech¬ 
niques facilitated development of transfer-function literal approximate 
factors for the F-4 in high AOA, non-zero-sideslip flight conditions. 

C. PROPOSED KEL-F-8785C MODIFICATIONS 

Based upon results of this program it is recommended that modifica¬ 
tions to MIL-F-8785C requirements be considered in three areas. 

I. Departure free Controlled Flight 

The current paragraph (3.4.2.2.1) simply requires that "all Classes 
of aircraft be extremely resistant to departure from controlled flight, 
post-stall gyrations and spins....The aircraft shall exhibit no uncom¬ 
manded motion which cannot be arrested promptly by simple application of 
pilot control." 

The "uncomoanded motio i“ requirement should be strengthened to 
require that, for roll control input within the service flight envelope 
(SFE) the following stability axis parameter values shall be no more 
negative than: 


Unaugmented airframes LCDP 



- 0.01 


Augmented airframe: 1/T^ > -0.5 

An alternate requirement might be that the aircraft should exhibit 
no aperiodic uncomoanded motion which exceeds 20 deg/sec and cannot be 
arrested promptly by simple application of pilot control. [The value 
selected here is based on a rovtgh average of the simulation pilots' 
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commentary as to their definitions of departure, the spiral divergence 
mode limit of Reference 40, and the 1/T^ limit above.] 

Since the SFE specifically excludes Btall and departure, the above 
criteria would automatically cover high AOA, pre-stall flight. Other 
potential locations for such wording are paragraphs 3.4.2.1.1, Stall 
Approach, and 3.4.2.1.2, Stall Characteristics. 

2. Stall Definition 

The current paragraph (3.4.2.1) allows definition of stall to be 
based on abrupt uncontrollable pitching, rolling, or yawing; or 

intolerable buffet. The results of our piloted simulation Indicate that 
any abrupt aperiodic rolling or yawing motion which octurB without being 
preceded by noticeable "g-break" is considered to be a departure, not 
a stall, and results in severely downgraded flying qualities. Thus, 
abrupt roll or yaw motion should be deleted as an allowable definition 
of stall and should not occur prior to stall. If such characteristics 
cannot be achieved with the airframe alone, then the flight control 
system should prevent reaching the AOA at which the abrupt rolling or 
yawing motion is obtained. 

3. Departure Warning 

The current paragraphs (3.4.1.1, 3.4.2.1.1.1, and 3.4.2.1.1.2) 

covering warning pertain mainly to stall. But, since abrupt uncommanded 
roll or yaw is also considered to define stall, the Paras. 3.4.2.1.1.1 
and 3.4.2.1.1.2 requirements, in effect, define departure warning. 
These requirements are based upon fixed margins of V Btall or CL 8tflll and 
do not take into account the severity of any post-stall/departure mo¬ 
tion, uncontrolled altitude loss, or mission phase. 

Results of this simulation indicated that stall/departure warning 
margin should be related to the severity of any uncommanded motion or 
PSG. Time and altitude loss prior to recovery were heavily dependent 
upon the initial departure severity; pilot commentary indicated the 
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overall departure ratings obtained were heavily influenced by alti¬ 
tude loss and mission phase. Quite obviously, consequences of stall/ 
departure In a ground attack phase can be more catastrophic than In 
high-altitude air combat. Unfortunately, insufClcier.'; data ware ob¬ 
tained to establish a specific departure warning crlt ^/lon. 

4. Compliance Demonstration 

The current specification (Para. 4.1) allows compliance with all 
requirements of Section 3 to be demonstrated through analyele. Although 
compliance with many of the requirements will be demonstrated by simula¬ 
tion, test, or both, this analysis and simulation program has demon¬ 
strated that linearized, frozen-point analytic results may be' aidleading 
or extemely difficult to interpret for flight regions where the airframe 
may exhibit large static aerodynamic cross-coupling momenta. Therefore, 
It is recommended that the requirement be changed so that compliance 
with paragraphs pertaining to stall/departure must be demonstrated via 
simulation and/or flight. 
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